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It has been shown that different Abelian and non-Abelian fraction quantum Hall states can be 
characterized by patterns of zeros described by sequences of integers {S^}. In this paper, we will 
show how to use the data {S a } to calculate various topological properties of the corresponding 
fraction quantum Hall state, such as the number of possible quasiparticle types and their quantum 
numbers, as well as the actions of the quasiparticle tunneling and modular transformations on the 
degenerate ground states on torus. 
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I. INTRODUCTION 

Characterization and classification of wave functions 
with infinite variables are the key to gain a deeper un- 
derstanding of quantum phases and quantum phase tran- 
sitions. For a long time, people have used symmetry 
properties of the wave functions and the change of those 
symmetries to understand quantum phases and quantum 
phase transitions. However, the discovery of fractional 
quantum Hall (FQH) states 1,2 suggests that the symme- 
try characterization is not enough to describe different 
FQH states. Thus FQH states contain a new kind of 
order - topological order. 3,4 Topological order is so new 
that the conventional mathematical tools and the lan- 
guage developed for symmetry breaking orders are not 
adequate to describe it. Thus, looking for new mathe- 
matical framework to describe the topological order be- 
comes an important task in advancing the theory of topo- 
logical order. 5-8 Intuitively, topological order correspond 
to pattern of long-range quantum entanglement in the 
ground state, as revealed by the topological entanglement 
entropy. 9,10 Characterizing topological order is character- 
izing patterns of long-range entanglement. 

In a recent paper, 11 the pattern of zeros is introduced 
to characterize and classify symmetric polynomials of in- 
finite variables. Applying this result to FQH states, we 
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find that the pattern of zeros characterizes and classifies 
many Abelian and non-Abelian FQH states. In other 
words, the pattern of zeros characterizes the topological 
order in FQH states. 4 This may lead to a deeper under- 
standing of topological order in FQH states. 

More concretely, a pattern of zeros that characterizes 
FQH ground state wave functions is described a sequence 
of integers {S a }. Describing pattern of zeros through a 
sequence of integers {S a } generalizes the pseudo poten- 
tial construction of ideal Hamiltonian and the associated 
zero-energy ground state. 12-17 All the topological prop- 
erties of FQH states should be determined by the data 
{S a }- In this paper, we will show how to use the data 
{S a } to calculate those topological properties. They in- 
clude the quantum number of possible quasiparticles and 
the ground state degeneracy on torus. 18,19 We also study 
the actions of quasiparticlc tunneling and the modular 
transformation on the degenerate ground states. 20 The 
main results of the paper are summarized in Section V. 

II. FUSION AND PATTERN OF ZEROS 

In Rcf. 11, we have introduced the pattern of zeros of 
a wave function by bringing a variables together. In this 
section, we will review the pattern of zeros within such an 
approach, but from a slightly different angle. For readers 
who are interested in the new results can go directly to 
section III. 



A. Fusion of a variables 

To bring a variables in a symmetric polynomial 
$({zi}) together, let us view zi,--- ,z a as variables 
and z a+ i, z a+ 2, • • • as fixed parameters. We write 
$(zi,z 2 ,- • •) as 

$({zj) = $(zi,--- ,z a ;z a+1 ,z a+ 2,---) 

where $(zi, • • • ,z a ; z a +i, z 0+2 , • • • ) is a symmetric poly- 
nomial of zi, • • • , z a parameterized by (z a+ i, z a+2 , ■••). 
Now let us rewrite 

z i = A£ + zW z^ = Zl + '" + Za , 

a 

a 

^& = 0, i = !,■■■, a, 

»=i 

and expand $(A£i + z (a \ ■ ■ • , A£ Q + z^\ z a+1 , z a+2 , • • • ) 
in powers of A: 

$(Aa + *<«*), ■■-,*& + *<°>; z a+1 , z a+2 , • • • ) 

OO 

= £ A m P m [*(°\ (&,■■■ ,£ a );z a+1 ,z a+2 ,...] 

m=0 

where P m [z^ a) , (£i, ■ ■ ■ ,£a); z a +i, z a+2 , ■••] is a homoge- 
nous symmetric polynomial of , £ a of degree m. 



P m [z (a) , (£i, • • • ,£ a );z a+ i,z a+2 , • • •] is also a polynomial 
of z( a ) and z a +i, z a+2 , • • • . The minimal power of z^ is 
zero. 

The polynomial P m [z {a \ • • • , Z a +l, Z a +2, ' ' ' \ 
is derived from <E>({zi}) after fusing a z^-variables 
into a single z^ variable. Different choices of 
£i , ■ • • , £ a represent different ways to fuse z\ , ■ ■ ■ , z a into 
z( a \ As a polynomial of z^ a \ z a +\, z a + 2 , • • • , differ- 
ent ways of fusion can produce linearly independent 

P m [z( Q ),(a,--- ,ia);Z a+1 ,Z a+2 ,---\. Thus ■■■,&) 

can be viewed as a label for those linearly independent 
polynomials of z( a \ z a+1 , z a+2 , ■■■ ■ 

P m [z {a \ • • • ,£ n ); z a+1 ,z a+2 , ■■■} is also a symmet- 
ric polynomial of £i, ■ ■ ■ , £ . Different z^ a \ z a+ i, z a+2 , • • • 
also lead to different symmetric polynomials of 
, £ a - Let if a ,m be the number of lin- 
early independent symmetric polynomials of £i, • • • ,£ a , 
P m [z (a) , (£i, • • • ,£ a );z a+ i,z a+2 , ■• •], produced by all the 
choices of z (a) , z a+1 , z a+2 , ■ ■■ . Let i<^™j(£i, ■ • • ,£„) be a 
basis of those symmetric polynomials of £i , • • • , £ a , where 
a ( a ) = 1,2,--- ,K a . m label those linearly independent 
symmetric polynomials of £i, • • • , £ a . Thus the different 
polynomials of z( a \ z a+1 , z a+2 , ■ ■ ■ labeled by (£i, • • • ,&,) 
can be written as 

P m [z (a \ (£l, ' ' ' ,£a); 2a+l,^a+2j ' ' '] 

a (a,m) = l 

We see that as different polynomials of 
z( a \ z 0+ i, z Q+2 , • • • labeled by £i,"->£a, 
Pm[z (a) , (6, ' ' ' ,^a);z a+ i,z a+2 , • • •] can also be la- 
beled by a( a > m \ F^J itn) (t;i, ■ ■ ■ , £ a ) can be viewed as a 
conversion factor between the two different labels, o^"'" 1 * 1 
and 

Let S a be the minimal total power of Zi , • • • ,z a in 
$({zj}). We see that K a ^ m = if m < S a , and 

if a EE K a , Sa * 0. 

Since the symmetric polynomial F^ m) (£i, • • • , £ a ) of or- 
der m — 1 must have a form 

^V(fi.-,6.)«E6 

which vanishes due to the ^ £j = condition, we find 
that 

Sa * 1. (1) 

Here we would like to introduce a unique-fusion con- 
dition. If i^ a ,m a = 1 for all a and allowed m a , then we 
say that the symmetric polynomial <j>({zj}) satisfies the 
unique-fusion condition. In this paper, we will mainly 
concentrate on symmetric polynomials <f>({zi}) that sat- 
isfies the unique fusion condition. In this case we can 
drop the a[ a ^ labels. 
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Note that Si is the minimum power of Zj in $({zj}). If 
S\ > 0, we can construct another symmetric polynomial 
with less total power of Zj 



$({z i })=$({z i }vn^ 



.Si 



Since Y\ i z i 1 is non-zero except at z = 0, $({2^}) and 
$({zj})/ 1 have the same pattern of zeros. Since 
the minimal power of Zj in <£>({zi}) is zero, without loss 
of generality, we will assume 



Si = 0. 



(2) 



For any symmetric polynomial, we obtain a sequence 
of integers Si, S2, • • • . Different types of symmetric 
polynomials may lead to different sequences of integers 
Si , S2 , • • • • Thus we can use such sequences to character- 
ize different symmetric polynomials. Since S a describes 
how fast a symmetric polynomial $({zi}) approaches 
to zero as we bring a variables together, we will call 
5i, 5*2, • • • a pattern of zeros. 

For the v = 1/q Laughlin state (3) 



i<3 



we find that S a is 



S a = 



qa(a — 1) 



(3) 



(4) 



For the Laughlin state, K a ,m a = or 1. Thus the Laugh- 
lin state is a state that satisfies the unique-fusion condi- 
tion. 



This allows us to use S a to calculate D a b- It also implies 
that D a b does not depend on af\ and a^ a+b \ 
To derive (6), let us write $({zi}) as 

$(zi, • • • , Z a+b ; Za+b+l, Za+b+2, ' ' ' ) 

= x s a ~ x s b J- P(z^,a^;z^,a^;z a+b+1 ,---)x 
F a {a) (£1 , • • • , £ a )F a (.b) (£ a +i , • • • , S,a+b) H • (7) 

where 

z (a) = («! + ■■■ + z a )/a, z< 6) = (z a+ i + • • • + z a+b )/b, 
A& = z i — z< " a \ for i = 1, • • • , a, 
Mi — z i ~ z< " b \ for i = o + 1, ■ ■ ■ , a + b. 

In (7), • • • represents the higher order terms in A and A. 
Since in <&({zj}), the minimal total power of Zi, • • • ,z a 
is S a and the minimal total power of z a +i, ■ ■ ■ , z a+b is 
S b , as a result F aM (£1, • • • , £ a ) and F a(b) (£ a+1 , • • • , £ a+6 ) 
are homogenous polynomial of degree S a and S& respec- 
tively. Since the minimal total power of z\ , • • • , z a+ b is 

S a +b, thus the minimal total power of z^ and z^ b ' in 
P(z( a \ a(») ; z( fe ) , a< 6 > ; z a+6+1 , z a+6+2 , • • • ) is S a+6 - S a - 
Sa. This proves (6). 

(6) expresses D a b in terms of S a . We can also express 
S a in terms of D a b- Since Si = 0, we find that 



Sa+l — S a + -D a ,l 



which leads to 



B. Derived polynomials 



In the last section, we have obtained a new polynomial 
, • • • ) from <i>({zj}) by fusing a variables 
together. We can fuse many clusters of variables together 
and obtain, in a similar way, a derived polynomial 



P{{z. 



}) 



where zf^ is a type-a variable obtained by fusing a Zj- 
variables together. Since different ways to fuse z,'s into a 
z^ a ' may lead to different derived polynomials, the index 
cxf^ is introduced to label those different derived poly- 
nomials. 



As we let 



» 



z 



0>) 



, we obtain 



(«) .,( Q ).J b ) „W. 



P(zJ 
- 



Q: 



■)|,(. 



5 W= z (a+I.) ( 5 ) 



D a h is another set of data that characterize the symmet- 
ric polynomial $({zj}). The two sets of data, S a and -D b 
are related by 



D a b — S a +b — S a — Sb- 



(6) 



s a = j2 D b,i- 



(8) 



6=1 



n-cluster condition 



In Rcf. 11, we have introduced an n-cluster condition 
on symmetric polynomials $({zi}) to make polynomials 
of infinite variables to behave more like polynomials of 
n variables. Here we will introduce the n-cluster condi- 
tion in a slightly different way. A symmetric polynomial 
<3>({zj}) satisfies the n-cluster condition if and only if 

(a) the fusion of kn variables is unique: Kk n ,s kn = L 
where k =integer, 

(b) as a function of a /cn-cluster variable z^ kn \ the de- 
rived polynomial P{z^ kn \ z[ a \ • • • ) is non-zero if z( feTl ) ^ 

z- a ^ (ie no zeros off the variables z- a ^). 
Let 

m = D n) i, 
the condition (b) allows us to show that 
DknM = kam. 



(9) 
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From £>„,a = S a+n - S a - S n , we find that S a+n = S a + 
S n + ma which implies that 

k(k — 1) 

Sa+kn = S a + kS„ + TO71 h fcma. (10) 

(10) is the defining relation that defines the n-cluster 
condition on the pattern of zeros {S a }. We see that 
{S2, • • • , S n } determine the whole sequence {S a }- 



D. Occupation description of pattern of zeros 

There are other ways to encode the data {S a } which 
can be very useful. In the following, we will introduce 
occupation description of the pattern of zeros {S^}. We 
first assume that there are many orbitals labeled by 
I = 0, 1, 2, • • • . Some orbitals are occupied by particles 
and others are not occupied. The pattern of zeros {S a } 
corresponds to a particular pattern of occupation. To 
find the corresponding occupation pattern, we introduce 

l a = S a -S a -l, 0=1,2,---. (11) 

and interpret l a as the label of the orbital that is occupied 
by the a th particle. Thus the sequence of integers l a 
describe a pattern of occupation. l a is a monotonically 
increasing sequence: 

L+2 - la+1 = S a+ 2 + S a — 2S a +l > 

as implied by the condition (see Ref. 11) 

S a +b+c — S a +b — S a +c — Sb+c + S a + Sb + S c > (12) 

with b = c = 1. The n-clustcr condition (10) becomes 

L+n = l a + m. (13) 

So {li, • • • , l n } determine the whole sequence {l a }- Two 
sequences, ■ ■ ■ ,l n } and {5*2, • • • , S n }, have a one-to- 
one correspondence. Thus {S a } and {l a } arc faithful rep- 
resentations of each other. 

The occupation distribution can also be described by 
another sequence of integers n\. Here rii is the number 
of la's whose value is I. ni is the number of particles oc- 
cupying the orbital I. Both {/ a } and {n;} are equivalent 
occupation descriptions of the pattern of zeros {S a }. 

The occupation distribution {m} (or equivalently {l a }) 
that describes the pattern of zeros in <f> has a very physi- 
cal meaning. If we regard ni as the occupation number on 
an orbital described by one-body wave function </>; = z l , 
then the occupation distribution {n{\ will correspond to 
a many-boson state described by a symmetric polynomial 
${„ ( }. The two many-boson states $ and 3>{ n; } will have 
similar density distribution, in particular in the region 
far from z = 0. 



E. Ideal Hamiltonian and zero energy ground state 

For an electron system on a sphere with flux 
quanta, each electron carries an orbital angular momen- 
tum J = N<p/2 if the electrons are in the first Landau 
level. 12 For a cluster of a electrons, the maximum allowed 
angular momentum is aJ. However, for the wave func- 
tion $({zi}) described by the pattern of zeros {S a }, the 
maximum allowed angular momentum for an a-electron 
cluster in $({zi}) is only J a = a J — S a . The pat- 
tern of zeros forbids the appearance of angular momenta 
aJ—S a + l,aJ—S a +2,--- , a J for any a-electron clusters 
in *({*}). 

Such a property of the wave function $({zj}) can help 
us to construct an ideal Hamiltonian such that $ is the 
exact zero-energy ground state. 12-17 The ideal Hamilto- 
nian has the following form 

%u=e E (14) 

a a-clectron clusters 

where is a projection operator that projects onto 
the subspace of a electrons with total angular momenta 
aJ—S+1,-- - , a J and Pf la is a projection operator into 
the space H a . 

What is the space TL a l Let us fix z a +i, z a+ 2, • • • , and 
view the wave function $(zi,- • • ,Zn) as a function of 
zi,--- ,z a . Such a wave function describes an a-electron 
state. If $(zi, • • • ,zn) is described by a pattern of zeros 
{S a }, the a-electron state defined above can have a non 
zero projection into the space H a ,s a , where H a ,s a is a 
space with a total angular mentum a J — S a . However, 
different positions of other electrons may lead to different 
projections of the a-electron states in the space H a ,s a - 
H a is then the subspace of H a ,s a that is spanned by 
those states. ft a is the subspace of Ti. a ,s a formed by 
vectors that are orthogonal to 7i a . 

We see that by construction, <I> described by a pattern 
of zeros {S a } is a zero energy ground state of the ideal 
Hamiltonian -ff{s a }- However, the zero energy states of 
H{s a } are not unique. In fact, any state $ with a pattern 
of zeros {S a } will be a zero-energy eigenstate of -ff{s a } 
if S a > S a - (If S a = S a , we will further require that 

H a C H a .) 

But on sphere, H^' a } ma y have a unique ground state. 
If the number of electrons is a multiple of n: N — nN c , 
and the number of the flux quanta N$ satisfies 

2 9 

2J = N 4) = —^+m(N c -l), (15) 

then the state $ can be put on a sphere and corresponds 
to a unique uniform state with zero total angular momen- 
tum. Such a state may be the unique zero-energy ground 
state of the ideal Hamiltonian i?{s a } on the sphere. 

If we increase (but fix the number of particles N) , 
then -ff{5 a } will have more zero-energy states. Those 
zero-energy states can be viewed as formed by a few 
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quasiparticle excitations. Those quasiparticle excitations 
may carry fraction charges and fractional statistics. In 
the next section, we will study those zero-energy quasi- 
particles through the pattern of zeros. 



III. QUASIPARTICLES 

A. Quasiparticles and patterns of zeros 

If we let z\ , ■ ■ ■ , z a approach in a ground state wave 
function $({zi}), we obtain 

*({*})U-o = A 5 »P(6, z„+i, •••) + ••• 

where Z{ = A£j, i = 1, • • • , a. The pattern of zeros {S a } 
indicates that the z = point is the same as any other 
point, since we will get the same pattern of zeros if we let 
Zi, ■ ■ ■ , z a to approach to any other point. Thus the pat- 
tern of zeros {S a } describes a state with no quasiparticle 
at z = 0. 

If a new symmetric polynomial <& 7 ({zi}) has a quasi- 
particle at z = 0, $ 7 ({zi}) will have a different pattern 
of zeros {S TtU } as z\, ■ ■ ■ ,z a approach 0: 



$ 7 ({zi}) - A^»P 7 (a,--- ,£ a ;z a+ i, •••) + ■ 



(16) 



In other words, the minimal total powers of z\, ■ ■ ■ , z a in 
<j> 7 ({zi}) is S-y-a- Thus we can use {5 7;a } to characterize 
different quasiparticles. Here the index 7 labels different 
type of quasiparticles. We would like to stress that the 
quasiparticle wave function <& 7 ({zi}) is still described by 
the pattern of zeros {S a } if we let z\, ■ ■ ■ , z a approach to 
any other point away from z — 0. 

Just as {S a }, {S-y-a} should also satisfy certain condi- 
tions. Here, we will consider quasiparticle states <& 7 ({zi}) 
that have a zero energy for the ideal Hamiltonian H^g a y. 
The zero-energy condition on quasiparticles requires that 
S T>a should satisfy 



S~/:a ^ S a . 



(17) 



Although both the ground state $ and the quasiparticle 
states $ 7 have a zero energy, the ground state has the 
minimal power of Zi while the quasiparticle states have 
higher total powers of Zj. 

In the following, we will discuss other conditions on 



S T , a . Let P 7 ({ 



\ }) be the derived polynomial ob- 



tained from $ 7 ({z i }). First we let z( a ) -» m P~,: 

P- ( (z^\a^;z^ ia ^;---) 
~ {z^) D ^P', ( {a {a ^,z^\ a ^ ;•••)• 

Then we let z^ — > and let P 7;a ,o be the order of zeros 
in (a( Q ) ; z^ , aW ; • • • ) as z^ 0: 

J*(a<°>; *< 6 >,a< 6 >; ■ ■ ■ ) ~ ( z ( i ))^M P ''( a W ;a W ; . . . ). 



Y;fl • 

n x 



FIG. 1: As a function of z (c) , P 7 (a (a) ; z (6) , a (ft) ; z (c) , a (c) ; ■ • • ) 
has a _D 7;ajC th order zero as — > 0, and a -Dj,, c th order zero 
as — ► z' 6 '. The crosses mark the positions of other zeros 
that are not at z = and not at z* 6 - 1 . 



We have 

+ D T , a . b + S b , (18) 
which leads to a condition on S TiU 

Sj;a+b ~ S^- a — Sb > 0. 

To get more conditions on S* 7;a , let us consider fusing 
three variables z^ a \ z^ h \ and z^ to 0. We may first fuse 
z< a ) to to obtain P 7 (a^ ;z^, ; z< c ) , a ( c ) ;•'••)• Then 
we fuse z^ to to obtain P^ (a^ ; ; z^ , ; ■ ■ ■ ) . 
Last we fuse z^ to 0. In this case, we obtain a D 7:a +6 !C th 
order zero as z^ -> in ^'(aW; qW; z( c \ qW; • • • ). We 
also know that P 7 has a P , c th order zero as z( c ) — > z^ b ' 
and a P 7:a c th order zero as z( c ) — > (see 1). If z( b ) is 
very close to 0, we find that P 7;a+ b iC is the sum of £> 7;a>c , 
Pfc, c , and the zeros close to z = but not at z = and 
at z^ b \ This way, we find 



-^7;a+6,c ^ ^7;a,c ~t~ Po,< 

which implies that 



(19) 



^[a+fc+c — 5*7:0+0 — 5 7:a + c Sb+c + S~ f , a + 5*o + S c > 0. 

(20) 

(20) generalizes the condition (12) on the pattern of zeros 
of the ground state. 

The n-cluster condition implies that, as a function of 
z (n \ P 7 (z("),zf a) ,---) is non-zero if z^ ^ z\ a) and 
z(") 7^ 0. Therefore, the inequality (19) is saturated when 
c = n: 

which implies that (setting a = 0) 

S^fc+n = S-yjo + P 7; o,n + Sfc+ri — Sb 

= + 5„ + P 7; o,n + mb, 

or 



5 7;a+ fe„ = S Tia + k(S n + ma + D 1 , . n ) + ran 



k(k - 1) 



Since Sj — 0, we have S j;n = S n + D 7; o,n- Thus 

A;(A; — 1) 

= 5 7;a + k(S T . n + ma) + ran . (21) 

for a > 0. This is the n-cluster condition on S^- a . 



6 



B. Solutions of patterns of zeros for quasiparticles 

The ground state of the ideal Hamiltonian (14) is a 
state described by a pattern of zero {S'a}- The zero- 
energy quasiparticles above the ground state <& can be 
characterized by a pattern of zeros {S i;a }. If $ is a state 
with the n-cluster form, then all the 5 7;a 's are deter- 
mined from S-y-i, ■ ■ ■ , S-y-n (see (21)). Thus zero-energy 
quasiparticles arc labeled by n integers: S Ti i, ■ ■ ■ , S T . n . 
From the discussion in the last section, we find that those 
integers should also satisfy 

S T: a+b — Sj-a — Sb > 0, (22) 

S'y-a+b+c ~ 5*7:0+6 — S^y-a+c ~~ Sb+c + S~ r .a + Sb + S c > 0. 

Numerical experiments suggest that the solutions of (22) 
always satisfy (17). Thus the condition (17) is not in- 
cluded. 

For a given pattern of zeros {S a } describing a sym- 
metric polynomial <& of n-cluster from, (22) has many 
solutions that satisfy the n-cluster condition (21). Those 
solutions can be grouped into equivalent classes. 

To describe such equivalent classes, we need to use the 
occupation description of the pattern of zeros {S T>a }. We 
introduce 

lj :a Sj- a Sj- a —i, a — 1, 2, • . 

and interpreted Z 7;a as the label of the orbital that is 
occupied by the a th particle. Thus the sequence of in- 
tegers l 7 - a describes a pattern of occupation. Wc note 
that l~[ ]a +i — D~ ( ; a .i. (19) implies that / 7;a +i > l-y;a- The 
occupation distribution can also be described by another 
sequence of integers n 7; /. Here n 7; ; is the number of i 7;a 's 
whose value is I. n 7; ; is the number of particles occupy- 
ing the orbital I. Both {^ 7:a } and {n 7; /} are equivalent 
occupation descriptions of the pattern of zeros {S T . a }- 
The n-cluster condition on {S-y- a } implies that 

^7;a+n = ^7;a + W^i (^) 

for a > 1. This implies that occupation numbers {n 7; /} 
satisfy 

n 7; ;+m = rij-i (24) 

for a large enough I. Thus the occupation numbers n 7: ; 
become a periodic function of I with a period m in large 
/ limit. In that limit, each to orbitals contain n particles. 

Now consider two quasiparticles 71 and 72 described by 
two sequences of occupation numbers {n] 1 } and {nj 2 }. 
If n 7i; ; = n 72: ; for large enough I, then the density distri- 
butions of two many-boson states <E> 71 and <I> 72 only differ 
by an integral number of electrons near z = 0. Thus the 
two quasiparticles 71 and 72 differ only by an integral 
number of electrons. This motivates us to say that the 
two quasiparticles to be equivalent. 

Numerical experiments suggest that each equivalent 
class is represented by a simple occupation distribution 
n 7: ; that is periodic for all I > 0: 

n 7 :i +n = n T: i, l>0. (25) 



We call such a distribution the canonical distribution for 
the corresponding equivalent class. For a canonical oc- 
cupation distribution, there are n particles in every unit 
cell (I = 0, • • • , m — 1), (/ = m, ■ ■ ■ , 2m — 1), 

We can obtain all the equivalent classes of the quasi- 
particles by finding all the canonical distributions that 
satisfy (22). The equivalent classes of the quasiparticles 
correspond to fractionalized excitations. 



IV. TOPOLOGICAL PROPERTIES FROM 
PATTERN OF ZEROS 

A FQH state characterized by a pattern of zeros {S a } 
can have many topological properties, such as quasipar- 
ticle quantum numbers, 2,21 ground state degeneracy on 
compact space, 3,19,22 and edge excitations. 23 In this sec- 
tion, we are going to calculate some of those topological 
properties from the data {S a }. 

A. Charge of quasiparticles 

We have seen that a quasiparticle excitation labeled 
by 7 are characterized by a sequence of integers: {5 7;o }. 
We would like to calculate the quantum numbers of the 
quasiparticle from {S-y- a }- 

To calculate the quasiparticle charge, we compare the 
occupations m that describe the pattern of zeros {S a } 
of the ground state $ and the occupations n T j that de- 
scribe the pattern of zeros {S Tta } of a quasiparticle state 
<I> 7 . We divide I = 0, 1, 2, • • • into unit cells each contain- 
ing to orbitals: I = 0, • • • ,m — l;m, ■ ■ ■ , 2to — 1; ■ ■ ■ . n/ 
and n 7; ; contain the same number of particles in the kth 
unit cell if k is large enough. Since n; is a distribution 
that corresponds to zero quasiparticle charge, we might 
think that the quasiparticle charge corresponding to the 
distribution n 7: ; is given by 

km 

1=0 

in large k limit. However, this result is incorrect. Al- 
though n; and n 7; / contain the same n particles in the 
kth unit cell (for a large k), the "centers of the mass" of 
the two distributions in the fcth cell are different. The 
shift of the "centers of mass" is given by 

km—l 

Shifting the "center of mass" by to is equivalent to 
adding/removing n particles. Thus the total quasipar- 
ticle charge is given by 

km _^ km—l 

Qi = ^2( n i - n i;i) - — ( n l- n r,l)h (26) 

1=0 l—km—m 
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FIG. 2: The 2 J + 1 orbital on a sphere form a angular mo- 
mentum J representation of the 0(3) rotation. 



or more generally the Berry's phase of the quasiparticle 
is 

el=nJ2(n T> i-m)(l-J). (29) 

where the upper bound of the summation X^=o ^ s 
roughly at I = = 2J. From the part of B that is 
linear in = 2 J, we can recover (26) for the charge of 
the quasiparticle. 



for a large enough k. Note that in the above definition, a 
charge +1 corresponds to an absence of an electron. For 
a canonical occupation distribution satisfying (25), the 

first term Y^i=o( n i ~ n id) vanishes. 

Since the two descriptions of occupation distributions, 
{^7;a} and {n 7; i}, have a one-to-one correspondence, we 
can also express Q 7 in terms of {/ 



according to (23), {I 
sequence {£ 7;a }): 



7;ii • • • i l~/;n} (note that, 
7;i> ■ ■ " ! ^7;n} determines the whole 



1 " 



(27) 



There is another way to calculate the charge of a quasi- 
particle. We can put the quasiparticle state $ 7 on a 
sphere with flux quanta. If we move the quasiparticle 
around a loop that spans a solid angle fi, the quasiparti- 
cle state $ 7 will generate a Berry's phase 



Q 7 + 0(N°) 



(28) 



The part in 0b that is proportional to allows us to 
determine the charge Q 7 . 

From the Berry's phase 0b(Q) (as a function of the 
solid angle ft), we can find out the angular momentum 
J q of the quasiparticle 

j q = B (n)/n, 

(note that 0b oc f2.) The occupation distribution ni for 
the ground state $ describes a trivial quasiparticle with 
zero charge and zero angular momentum. The occupa- 
tion distribution n 7; ; for the quasiparticle state <I> 7 de- 
scribes a non-trivial quasiparticle with a non-zero angu- 
lar momentum J q . Since the orbital 0/, when put on a 
sphere with flux quanta, is identified with an angular 
momentum eigenstate (J, J z ) — (-^-,1 — J), the quasipar- 
ticle state $ 7 is an eigenstate of total J' ot (see Fig. 2). 
Since $ 7 is the state with maximal J* ot , the total angular 



momentum of the quasiparticle is J 



J tot — T 
y — » a 



Since the occupation distribution ni describes a state 
with J* ot = 0, the angular momentum of the quasiparti- 
cle is 



Jg = ^2( n ~/:i - n i)( l - J ), 



B. Orbital spin of quasiparticles 

We have seen that the Berry's phase of the quasipar- 
ticle contains a term linear in = 2 J that is related to 
the quasiparticle charge. The Berry's phase also contains 
a constant term which by definition determines the or- 
bital spin S° sp of the quasiparticle. 4 ' 24,25 More precisely, 
in the large = 2 J limit, we have 



1 

n 



= JQ 7 + S° sp + O^J- 1 ). 



(30) 



Therefore, to calculate S° sp , we need to evaluate (29) 
carefully. (29) is not a well defined expression since the 
upper bound of the summation X)z=o i s n °t given pre- 
cisely (on purpose). So to evaluate (29), we first regulate 
(29) as 



6 B = ClJ2(n Til - m)(l - J)e 



a I 



(31) 



and then take the small a and large J limit with a ~ 
1/J. The key difference between (29) and (31) is that 
(31) has a soft (or smooth) cut-off near the upper bound 
of the summation J2i=o- After evaluating the regulated 
summation (31) in appendix 1, we find that the orbital 
spin of the quasiparticle is given by 



2=0 



2m J 



(32) 



(32) is valid only for canonical occupation distributions 
that satisfy (25). In terms of l 7ia , we can rewrite (32) as 



S° sp 



7;a 



i 2 

7;a 



ml a + 1\ 



2m 



(33) 



However, (33) is more general. It applies to generic quasi- 
particles, even those whose occupation distributions do 
not satisfy (25). 

We would like to stress that cqn. (33) is obtained 
through an untested method. Its validity is confirmed 
only for Abclian quasiparticles (see appendix 2). Inde- 
pendent confirmation is needed for more general cases. 
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C. Ground state degeneracy on torus 

A FQH state has a topological ground state de- 
generacy on torus, which is robust against any local 
perturbations. 3,19 Such a topological ground state de- 
generacy is part of the defining properties of topological 
orders. 

According to topological quantum field theory, 26-28 the 
topological ground state degeneracy is equal to the num- 
ber of quasiparticle types. To be precise, two quasiparti- 
cles are regarded equivalent if they differ by a multiple of 
electrons. Thus, the topological ground state degeneracy 
is equal to the number of equivalent classes of quasipar- 
ticles introduced in section IIIB. 

To understand such a result, let us consider a quasipar- 
ticle state $ 7 on a sphere (which is a zero-energy state 
of the ideal Hamiltonian) . We stretch the sphere into a 
thin long tube. The state $ 7 remains to be a zero-energy 
state in such a limit. According to Ref. 29-31, the FQH 
state $ 7 in such a limit becomes a charge-density-wave 
(CDW) state characterized by a particle occupation dis- 
tribution among the orbitals on the thin tube. We ex- 
pect such an occupation distribution is given by n 7; ; in 
the large I limit. We see that n 7: ; in the large I limit is 
a CDW state in the thin cylinder limit that is compati- 
ble with the zero-energy requirement. This suggests that 
the quasiparticle types (determined by n 7 -i in the large I 
limit) and the zero-energy CDW states on thin cylinder 
(given by n 7; ; in the large I limit) are closely related. 

The above result allows us to construct zero energy 
ground states on a torus. Let us consider a torus with 
N<p = N c m flux quanta. There are orbitals on such a 
torus. Those orbitals are labeled by I = 0, 1, 2, • • • , — 
1. Now let us consider aJV = A c n-electron FQH state on 
such a torus. The FQH state is described by a pattern 
of zeros {S a } of n-cluster form and has a filling fraction 
v = n/m. What are the degenerate ground states of such 
a FQH state on the torus? 

According to Ref. 29-31, the zero-energy ground state 
wave functions of a FQH state in the thin cylinder limit 
can be described by certain occupation distribution pat- 
terns (or certain CDW states). The above discussion on 
sphere suggests that the canonical distributions n 7; ; for 
the quasiparticles are just those occupation distributions. 
(Note that the canonical distributions fill each m orbitals 
with n electrons and give rise to very uniform distribu- 
tions.) Thus each canonical distribution n 7; ; gives rise 
to a TV-electron ground state wave function on the thin 
torus which corresponds to a degenerate ground state. 
We see that the canonical occupation distributions n 7; ; 
characterize both the degenerate ground states on torus 
and different types of quasiparticles. This explains why 
the ground state degeneracy on torus is equal to quasi- 
particle types. We also see that we can use m integers 
n 7; ;, I = 0, 1, • • • , m — 1, to label different degenerate 
ground states. 

The condition on the CDW distributions, n 7; ;, that 
correspond to the zero energy ground states on thin torus 



can be stated in the translation invariant way. Using 
^ 7; a's, we can rewrite the second expression of eqn. (22) 
as 

c 

^(^ 7 ;a+6+fc - l-f.a+k) > Sb+c - Sb - S c = D h .c- (34) 
k=l 

To understand the meaning of eqn. (34), let us set c = 1 
in the in eqn. (34). In this case, eqn. (34) requires that a 
zero energy CDW distribution must satisfy the following 
condition: any groups of b electrons must spread over 
Dbs + 1 orbitals or more. 

After knowing the one-to-one correspondence between 
the quasiparticle types and the degenerate ground states, 
we like to ask which quasiparticle type correspond to 
which ground states? Since both quasiparticle types and 
the degenerate ground states are labeled by the canonical 
occupation distribution, one may expect that a quasi- 
particle labeled by n 7; ; will correspond to the ground 
state labeled by n 7; /. However, this does not have to the 
case. In general, a quasiparticle labeled by n 7; / may cor- 
respond to the ground state labeled by n~ i; i + i s . Later, 
we see that the shift l s is indeed non-zero. So we will 
denote the ground state wave function that corresponds 
to a quasiparticle 7 as ${„ }({Xi}) or more briefly 
as$ 7 ({XJ). 

D. Quantum numbers of ground states 

The Hamiltonian of FQH state on torus has certain 
symmetries. The degenerate ground state on torus will 
form a representation of those symmetries. In this sec- 
tion, we will discuss some of those representations. 

1. The Hamiltonian on torus 

First, let us specify the Hamiltonian of the FQH system 
more carefully. The kinetic energy of the FQH system is 
determined by the following one-electron Hamiltonian on 
a torus (A 1 , A 2 ) ~ (A 1 + 1, A 2 ) ~ (A 1 , A 2 + 1) with a 
general mass matrix: 

Hx= -\ £ {^- iAi MAi- iA ^ (35) 

where g is the inverse-mass-matrix: 

g(r) = ( T " + % ~p\ (36) 

V T v T v ) 

and 

(A 1 ,A 2 ) = (-2irN 4> X 2 ,0) (37) 

gives rise to a uniform magnetic field with flux quanta 
going through the torus. The state in the first Landau 
level has a form 

<HA\A 2 ) = f(X 1 +TX 2 )e iirN * T ^' (38) 
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FIG. 3: The circular orbital wave function 4> {l) (X 1 , X 2 ) in the 
thin cylinder limit. A think line marks the positions where 
^Wp^X 2 ) is peaked. 



where r = t x + \r y and f{z) is a holomorphic function 
that satisfies the following periodic boundary condition 

f(z+a+bT) = /( z ) c -i"rfc 2 A^-i2,r&AV ^ b= integers. 

(39) 

The above holomorphic functions can be expanded by 
the following basis wave functions: 



/W(z|r)=^e lN * 



iw-(N tt ,k+l) 2 + i2TT(N 4 ,k+l)z 



(40) 



where I — 0, • • • , — 1. The corresponding wave func- 
tions 

^ l \x\X 2 \r) =/W(X 1 +rX 2 |r)e i7rJV ^ 2 ) 2 (41) 

are orbitals on the torus. 

The thin cylinder limit 29-32 is realized by taking r y — > 
oo. In such a limit 

^ l \X\X 2 \r) = ^(X 1 +rX 2 |r)e i ^( x2 ) 2 

_ c ^TN 4 ,(k+ 7 ^+X 2 f + i2-K(N 4 ,k+l)X 1 
k 

^ ^iivTN 4 ,[X 2 -(l~ 7 ^)] 2 +i2Tr(l-N <t ,)X 1 

We see that in the thin cylinder limit, the orbital wave 
function <j>^ (X 1 , X 2 \t) is non-zero only near X 2 = 1 — 
^ (see Fig. 3). 



2. Translation symmetry 

First let us consider the symmetry of the Hamilto- 
nian (35). Since the magnetic field is uniform, we expect 
the translation symmetry in both X 1 - and X 2 -dircctions. 
The Hamiltonian (35) does not depend on X 1 , thus 



■4 



l (d 1 ,0) HT (d 1> 0) 



= H, 



T {dl ,o) = e d ^ 



But the Hamiltonian (35) depends on X 2 and there seems 
no translation symmetry in X 2 -dircction. However, we 
do have a translation symmetry in X 2 -direction once we 
include the gauge transformation. The Hamiltonian (35) 
is invariant under X 2 — > X 2 + d 2 transformation followed 
by a e^N^xi + ij, £/(i)_ gauge transformation: 



T (o,d 2 ) HT (o.d 2 ) - H , T (o,d 2 ) - e 1 



2irN 4 ,d 2 X 1 + i<p e d 2 d x 2 



In general 

TjHTd = H Td = G i27rN 't> d 2 xl + i ' !rN 4> d i- d 2 e <i-9 
where we have chosen the constant phase in Td 

4> = nN^dxdy (43) 

to simplify the later calculations. The operator Td is 
called magnetic translation operator. So the Hamiltonian 
(35) does have translation symmetry in any directions. 
But the (magnetic) translations in different directions do 
not commute 



T d Td> = c l27TN 't( dld '2- d 2 d 'i)Td>Td- 



(44) 



and momenta in X 1 - and X 2 -directions cannot be well 
defined at the same time. 

However, when is an integer, TVq) and T(o,i) do 
commute and the wave function <j>^ in the first Landau 
levels satisfies 

T (uo) ^ =T (0 ^ =<f>^. 



Therefore <j>^ is a wave function that lives on the torus 
(X\X 2 ) ~ (X 1 + l,X 2 ) ~ (X\X 2 + 1). 

On a torus, the allowed translations are discrete since 
those translation must commute with T( 10 ) and T(o,i). 
The smallest translation in X 1 and X 2 directions are 
given by 



(45) 



which satisfies 

TiT 2 = c i27r/Ar *T 2 7i. 
We also find that (see appendix 3) 

The many-body ground state wave functions 
$ {n7iI+Ia} ({Xi}) = * 7 ({Xi}). of the FQH state 
on torus are labeled by the canonical occupation distri- 
butions (see (25)). In the thin cylinder limit t v — > oo, 
the many-body ground state wave functions 1+I j 
become the CDW wave functions described by the oc- 
cupation distributions {« 7: ;+; 3 } where there are n 7; ; + ; s 
electrons occupying the orbital <f>W . This allows us to 
obtain how ${ n7 ., +! }'s transform under translation: 



Ti<I>r i = e i27r ^*=* ™Tii+i»V^(t). 



- {n T ;l+i 3 } 

ATi — 1 

— c i27r E;=o "TiiC-is)/^^ 
= c i2^S 1 "7;-(^+ ''*' 2 !: "" )$ 



I.) 
{™^;! + i s }- 



If we choose 



Ns-m + 
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(which is always an integer since — mN c , m = even, 
and l max — S n — S n -i = even), we find that 

Ti$ 7 = c i2,r ^™" 1 ("-' i! ~" ( )'' /m $ 7 , = c i27rQ7 $ 7 , 
T 2 * 7 = $ Wi(1 _ 1)%m+I<} = $ 7 ,. (46) 

where we have used J2k=o ~ n ^max/2 (see (78)). 
Here 7' is the quasiparticle described by the canonical 
occupation distribution n^/.j — n 7 .a_u% m . We see that 
the eigenvalue of T\ is related to the charge of the corre- 
sponding quasiparticle Q 1 = Y^ih=o ( n i;l — n l)l/ m - The 
action of T 2 just shifts the occupation distribution by one 
step. The new distribution describes a new quasiparticle 
7'. 



3. Modular transformations 

The degenerate ground states on torus form a pro- 
jective representation of modular transformation 20 M = 

d) ^ &L(2, ^) ' wncre a,b,c,d E Z and ad — be = 1 

(see appendix 4) . The modular representation may con- 
tain information that completely characterize the topo- 
logical order in the corresponding FQH state. 
The modular transformations are generated by 

m t =(;;), Ms =(?" 1 )- (47) 

For every M € SL(2, Z), we have an invertible transfor- 
mation U(M) acting on the degenerate ground states on 
torus. U(M) 's satisfy 

U(M)U(M') ~ U(MM'). 

where ~ mean equal up to a total phase factor. The two 
generators Mr and Ms are represented as 

T = U(M T ), S = U(M S ). (48) 

M — — 1 is represented as 

C = U(-l) 

which is the quasiparticle conjugation operator (see (91)). 

T, S, T\ and T 2 have the following algebraic relation 
(see appendix 4) 

TTi = TiT, TT 2 = c i7T ™T 2 T 1 T, 

STi = T^, ST 2 = CTtf^S. (49) 

Since Mj = -1 and (M S M T ) 3 = -1, we also have 

S 2 = C, C 2 = 1, (ST) 3 = c i0 C. (50) 




FIG. 4: The tunneling processes A M and B (72) . 

E. Quasiparticle tunneling around a torus 

1. Quasiparticle tunneling operators 

To use the modular transformation to obtain the prop- 
erties of quasiparticles, it is useful to consider the fol- 
lowing tunneling process: (a) we first create a quasi- 
particle 71 and its anti quasiparticle 7^, then (b) move 
the quasiparticle around the torus to wrap the torus 
n\ times in the X 1 direction and n 2 times in the X 2 
direction, and last (c) we annihilate 71 and 7*. The 
quasiparticle-pair creation process in the step (a) is repre- 
sented by an operator that map no-quasiparticle-particle 
states to two-quasiparticle-particle states. The quasipar- 
ticle transport process in the step (b) is represented by 
an operator that map two-quasiparticle-particle states 
to two-quasiparticle-particle states. The quasiparticlc- 
pair annihilation process in the step (c) is represented 
by an operator that map two-quasiparticle-particle states 
to no-quasiparticle-particle states. The whole tunneling 
process induces a transformation between the degenerate 
ground states on the torus: 

i7)-< 7 :, ) „ 2) i7) = i7 / )(< 7 1 i ,i 2) )y, 

For Abclian FQH states, W^^ n ^ is always an invert- 
ible transformation. But for non-Abelian FQH states, 
w/ 71 ^ n may NOT be invertible. This is because when 
we create the quasiparticle-anti-quasiparticle pair in the 
step (a), the pair is in such a state that they fuse into 
the identity channel. But after wrapping the quasiparti- 
cle around the torus, the fusion channel may change and 
hence the pair may not be able to annihilate into ground 
in the step (c). In other words, the annihilate process 
in the step (c) represents a projection into the subspacc 
spanned by the degenerate ground states. 
Let (see Fig. 4) 

4(71) _ w (7l) - W {ll) 

A ~ "(i,o) > — (o,i) - 

A combination of two tunneling processes in the X 1 di- 
rection: A*- 71 ) and then A^ J2 \ induces a transformation 
^4(72)^4(71) on degenerate ground states. A combina- 
tion of the same two tunneling processes but with a differ- 
ent time order: A^ 2 ^ and then A^ J1 \ induces a transfor- 
mation A^^A^ on the degenerate ground states. We 
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FIG. 5: A general quasiparticle tunneling process. 



note that the two tunneling paths with different time or- 
ders can be deformed into each other smoothly. So they 
only differ by local perturbations. Due to the topological 
stability of the degenerate ground states 3 , local perturba- 
tions cannot change the degenerate ground states. There- 
fore A^ 71 ' and A^ 2 ^ commute, and similarly B^ 1 ^ and 
commute too. We see that A^'s can be simultane- 
ously diagonalized. Similarly, B^'s can also be simulta- 
neously diagonalized. Due to the 90° rotation symmetry, 
A^ and B^ have the same set of eigenvalues. But since 
A^ and B^ in general do not commute, we in general 
cannot simultaneously diagonalize A^ 7 ' and B^K 

The basis $ 7 = 3>{„ . !+! } described by the occupa- 
tion distribution n 7; ; + ; s on the orbitals <j>^ is a natural 
basis in which A 1 - 11 ) is diagonal. This is because the tun- 
neling process A^ 1 ^ docs not move quasiparticle in the 
X 2 direction, and hence does not modify the occupation 
distribution n 1 a + i 3 on the orbitals <f>^\ On the other 
hand, B^ does move quasiparticle in the X 2 direction 
and hence shifts the occupation distribution n 7; ;+; s that 
characterizes the ground states. Therefore B^ 1 ^ is not 
diagonal in the $ 7 basis. In particular, when acted on the 
state $o = ^{ni +i } that corresponds to the trivial quasi- 
particle, produces the state $ 71 that corresponds 
the quasiparticle 71 : 

B^)$o = 6 7l * 7l , (51) 

where 6 7l is a non-zero factor. (Note that 7 = 
0, 1, • • • , N q — 1 where 7 = corresponds to the triv- 
ial quasiparticle and N q is the number of quasiparticle 
types.) 



2. Tensor category structure in quasiparticle tunneling 
operators 

W/'jl n s is just a special kind of quasiparticle tunnel- 
ing. In general, we can create many pairs of quasipar- 
ticles, move them around each other, combine and/or 
split quasiparticles, and then annihilate all of them 
(see Fig. 4). In addition to the quasiparticle-pair 



creation process represented by a mapping from no- 
quasiparticle-particle states to two-quasiparticle-particle 
states and the quasiparticle-pair annihilation process rep- 
resented by a mapping from two-quasiparticle-particle 
states to no-quasiparticle-particle states, the more gen- 
eral tunneling process also contains the quasiparticle 
splitting process represented by a mapping from one- 
quasiparticle-particle states to two-quasiparticle-particle 
states and the quasiparticle fusion process represented 
by a mapping from two-quasiparticle-particle states to 
one-quasiparticle-particle states. We will use W{X) to 
represent the action of the whole tunneling process on 
the degenerate ground state where X represents the tun- 
neling path. W(^ n2 ) discussed above is just a special 
case of W(X). 

Due to the topological stability of the degenerate 
ground states W(A') can have some very nice alge- 
braic properties. However, in order for W(X) to 
have the nice algebraic structure we need to choose 
the operators that represent the quasiparticle-pair cre- 
ation/annihilation and quasiparticle splitting/fusion pro- 
cesses properly. We conjecture that after making those 
choices, W(X) can satisfy the following conditions: 



W{X local ) cx 1, (52) 




where i,j,--- = 0,l,---,N q — l label the N q quasiparticle 
types and i = corresponds to the trivial quasiparticle. 
Note that N q is also the ground state degeneracy on the 
torus. The shaded areas in eqn. (52) represent other parts 
of tunneling path. Note that there may be tunneling 
paths that connect disconnected shaded areas. 

Here X\ oca \ represents a tunneling path which has a 
compact support (ie there is no path in X\ oca \ that wraps 
around the torus). In this case W(X\ oca \) represents lo- 
cal perturbations that cannot mix different degenerate 
ground states on torus. 3 Thus W(X\ oca i) must be propor- 
tional to the identity. The second relation in eqn. (52) 
implies that the tunneling amplitude W(X) only depends 
on the topology of the tunneling path. A smooth defor- 
mation of tunneling path will not change W(X). 
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Strictly speaking, due to the quasiparticle charge 
and the external magnetic field, we only have 



W 



= c 10 W 



where 9 is 



a path dependent phase. However, we can restrict the 
tunneling paths to be on a properly designed grid, such 
as a grid formed by squares. We choose the grid such 
that each square contains k units of flux quanta where 
k is an integer that satisfies kQ y = integer for any 7. In 
this case, e 10 = 1 for the tunneling paths on the grid. 
Other relations in eqn. (52) arc motivated from tensor 
category theory. 6-8 

One may notice that the rules (52) are about planar 
graphs while the tunneling paths are three dimensional 
graphs. How can one apply rules of planar graphs to 
three dimensional graphs? Here we have picked a fixed 
direction of projection and projected the three dimen- 
sional tunneling paths to a two dimensional plane. 33 The 
rules (52) apply to such projected planar graphs. Also, 
the action of tunneling process can only be properly rep- 
resented by framed graphs in three dimensions, to take 
into account the phase factors generated by twisting the 
quasiparticles. Here we have assumed that there is a way 
to choose a canonical framing for each tunneling process, 
such that their projections satisfy eqn. (52). 

The coefficients (di, F^™ l ,cofj) must satisfy the follow- 
ing self consistent relations: 6,7 
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(53) 



fc=0 



Here i* is the antiquasiparticle of i, Vi is defined by Vi — 
Vi* = \fd~i, and 5 ijk is given by 



$ijk 



1, if i, j, k can fuse into a trivial quasiparticle, 
0, otherwise 



We would like to point out that here we only considered 
a quasiparticle fusion algebra ipiipj = ^2 k N^ipk that has 
a special property = 0, 1. 
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FIG. 6: (a): Two tunneling processes: A (7l) and A (72) . (b): 
The tunneling path of the above two tunneling processes can 
be deformed according to the fifth relation in eqn. (52) with 
i = 72, j = 72, k — 71, k = 71, m = 0, and n = 73. (c): 
The fifth, the fourth, and the third relations in eqn. (52) can 
reduce (b) to (c). 





x2 (a) 



FIG. 7: (a): The ground state $ 7 =o on a torus that corre- 
sponds to the trivial quasiparticle can be represented by an 
empty solid torus, (b): The other ground state $ 7 that cor- 
responds to a type 7 quasiparticle can be represented by an 
solid torus with a loop of type 7 in the center. 



3. Implications of tensor category structure 

As the first application of the above algebraic struc- 
ture, we find that (see Fig. 6) 

4(72) 4(71) - P 7 2 72 ° rl F 737 i* 7 2* 4(73) 
A A - 71 7i* 73 ° 7 3 ^72 71*0 A 

73 

~ U 71727 3 ^ 1 



(54) 



73 



We see that the algebra of forms a representation of 
fusion algebra ip 11 tp l2 = J2 l3 ^71 7273 ^73- The operators 
5(7) — SA^S^ 1 (see eqn. (69)) satisfy the same fusion 
algebra 



5(72)5(71) = Y^S-n-y^B^. 



(55) 



As the second application of the above algebraic struc- 
ture, we can represent the degenerate ground states on 
torus graphically. One of the degenerate ground state $0 
that corresponds to the trivial quasiparticle 7 = can 
be represented by an empty solid torus (see Fig. 7a). 
We denote such a state as |0) = <&o- Other degenerated 
ground states can be obtained by the action of the B^ 
operators 



\-r) = B^\o). 



(56) 
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FIG. 8: The graphic representation of ^7). 





FIG. 9: (a): The graphic representation of A^"> (b): The 
graphic representation of A 1 '"' 'I7). (c): The fifth and the 
fourth relations in eqn. (52) can deformed the graph in (b) to 
the graph in (c). The shaded area represents the hole of the 
torus. 



From eqn. (51), we see that $ 7 and I7) are related 



6 7 $ 7 



l7>- 



Since I7) is created by the tunneling operator \y) 
can be represented by adding a loop that corresponds to 
the operator to the center of the solid torus (see Fig. 
7b). 

1 7) is a natural basis for tensor category theory. The 
matrix elements of A^ and B^ have simple forms in 
such a basis. From eqn. (55), we see that 

B (7 2 ) jB (7 1 )|0) =B (7 2 )| 7l) = ^ 7i727 ,| 73) . 

73 

Therefore, in the basis I7), the matrix elements of B^ 2 ^ 
are given by the coefficients of fusion algebra 



R(72) _ X 
-°737l — "7l727 3 • 



(57) 



The action of A^ ) on I7) is represented by Fig. 8. From 
Fig. 9, we find that 



5 TC 

A (7')| 7) = _1V| 7>) 



(58) 



We see that A^ 1 ' is diagonal in the I7) basis. Let a 
be the eigenvalues of A^ 1 ) we see that 



(7') 
7 



,(7') 



cTC cTC 

7*7 ' _ l'7* 



cTC ■ 



(59) 



^7172 m e( l n - (58) i s tnc amplitude of two linked local 
loops (see Fig. 10). SZSL satisfies 



cTC _ cTC 

7172 7271 ' 



1 1 




FIG. 10: The amplitude of two linked local loops is a complex 
number S^^y 2 . 



Using the tensor category theory (52), one can also show 
that 7 

oTC cTC 

7l7 727 _ V- g ,5-TC = V^^TC B (7i) / 61 ) 

CTC 717273 737 773^7372 V U1 ^ 

°7 73 73 

which can be rewritten as 



oTC 
'-'07 



X , - V S* TC R(7i) (t^ ]- 1 ) , - A [ll) 
°77' — °773' D 7372VV' : ' ) J727' — /i 7 * 7 '*- 



In the operator form, the above becomes 

A (7i) = CS TC B^\S TC )- l C, 



(62) 



(63) 



where C is the charge conjugation operator C\"f} — I7*). 
We see that S TC can change the tunneling operator 5(71) 
to A^ 1 ' > . Eqn. (61) can also be rewritten as (assume S TC 
is invertible) 



E 



oTC cTC f7QTC\-l\ 
J 7l7°727" J > h 

cTC 

a 0-i 



= S, 



717273 



(64) 



We note that the above expression is invariant under 

cTC __. cTC f 

T\ and ±2 generate the translation symmetry of the 
torus. We expect that T\ and T 2 commute with the al- 
gebraic structure of the tensor category. Thus we expect 
that eqn. (46) keeps the same form in the new basis (7): 



ri|7> 



|7>, 



= IV). 



(65) 



where 7' is the quasiparticle described by the canonical 
occupation distribution n 7 >-i — « 7: (;-i)% m - 



4- Quasiparticle tunneling operators under modular 
transformation 



The transformations n n induced by quasiparti- 
(60) cle tunneling processes have certain algebraic relation 
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with the modular transformations U(M). From Fig. 
14, we see that the modular transformation Mj- changes 



s to 

(ni,n 2 ) VV (n 1 +n 2 ,n 2 )- 



TW, 



(71 ) 



w, 



(7i ) 



(111,112) (ni+7i2,n2) 



J. 



(66) 



Since the modular transformation M5 generates a 90° 
rotation, we find 



sw^ l] > = VF/ 7l) .s. 

(11,12) (-r»2,in) 



(67) 



Here Mg and Mg are given by (47) and T and S are given 
by (48). Also since the modular transformation M = — 1 
generates a 180° rotation, we find 



CW, (7l) , = VK/ 7l) ,C. 

(ni,n 2 ) (-nj ,-ni) 



(68) 



In terms of A^ 71 ^ and £?( 71 ) we can rewrite eqn. (66), 
cqn. (67), and eqn. (68) as 

j^(7i) _ ^4(71)^ 



SAM =B (7l) 5, 
55(71 ) c = CA M S,. 



(69) 



Let us assume that the set of the quasiparticle oper- 
ators can resolve all the degenerate ground states 
1 7), ie no two degenerate ground states share the com- 
mon set of eigenvalues for the operators Ah). In this 
case, the commutation relation TAh 1 ) = Ah^T implies 
that T is diagonal in the I7) basis. We will fix the over 
all phase factor of T by choosing Too = 1 • 

The operator C — S 2 is a charge conjugation operator. 
Its action on I7) is given by (see eqn. (92)) 

G\l) = |7*>- 

Compare A^) = S^B^S = CSB^S^C with 
(63), we find that S TC = SF where F is a diagonal 
matrix -F 77 ' = fyS^yi in the | 7 ) basis. Using 5 TC = SF, 
we can rewrite eqn. (59) and eqn. (64) as 



3717/7 _ ^717 _ (71) 
- So ~ ° 7 * 



Sojf- 



7 



'07 



As- 1 ), 



7 



'07 



-'717273 



(70) 



(71) 



We see that A^ and Bh) can be determined from 5*. 

Since S TC is symmetric, Sj^ = d 7 > and Sj^ = 1, 
once we know 5, we can use those conditions to fix F. 
Thus we can determine S TC from S. Once we know S TC , 
we can also calculate the CFT scaling dimension /i 7 for 
the quasiparticle 7 (see appendix 6) up to an integer: 7 



°7l72 ~~ °7l72 7 3 e 



i2-7r(/i T 



(72) 



The CFT scaling dimensions for the quasiparticles deter- 
mine the T matrix: 



f 77' 



i 2-77/1-, 



(73) 



F. Summary 



In this section, we calculated the charge and the or- 
bital spin of quasiparticle, as well as and the ground 
state degeneracy from the pattern of zeros {S a } of a 
FQH states. We also discussed the translation trans- 
formations, the modular transformations, and the trans- 
formations induced by the quasiparticle tunneling on the 
degenerate ground states. The algebra of those transfor- 
mation can help us to determine the quasiparticle statis- 
tics, quasiparticle quantum dimensions, and fusion al- 
gebra of the quasiparticles. In particular, we can use 
the algebra (49) and (50) to determine S, and then use 
eqn. (70), eqn. (71), and eqn. (72) to determine the quasi- 
particle tunneling operators, Ah) and Bh) , and quasi- 
particle scaling dimensions h 1 . The condition that the 
matrix elements of Bh) must be non- negative integers 
put further constraint on S. 



V. EXAMPLES 
A. Quasiparticles in FQH states 

In Rcf. 11, many FQH states are characterized and 
constructed through patterns of zeros. The pattern of 
zeros in FQH states can be characterized by a S'-vector 
5 = (to; 5 2 , • • • , S„), or a h-vector h = {^;hf, • • • , h s ^), 
or an occupation distribution n • • • n m _i). All those 
data contain information on two important integers n 
and to. n is the number of electrons in one cluster and 
to determines the filling fraction v = n/m. 

The FQH states constructed in Ref. 11 include the Z n 
parafermion states $-^ ; z„ introduced in Ref. 15,35. The 
patterns of zeros {S a 7 for those Z n parafermion states 
are obtained. The occupation distributions {n;} of those 
states agree with those obtained in Ref. 29,30,32. Ref. 11 
also obtained generalized Z n parafermion states <&„_ . z < fc ) 

and their patterns of zeros. $ n.. z « has a filling fraction 
v = n/m. Many other new FQH states and their patterns 
of zeros are also obtained in Ref. 11, such as the <5>^ ; c n 
and the <&™ ;Cri / Z2 states. 

Once we know the pattern of zeros of a FQH state, we 
can find all the quasiparticle excitations in such a state, 
by simply finding all S^ ia that satisfy (22) and (25) (note 
that n 7; ; in (25) are determined from S 7 - a ). From the 
pattern of zeros that characterizes a quasiparticles, we 
can find many quantum numbers of that quasiparticle. 
Here we will summarize those results by just listing the 
number of quasiparticle types in some FQH states. Then, 
we will give a more detailed discussion for few simple 
examples. 
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For the parafermion states $§ ; z"„ (m = 2), we find the numbers of quasiparticle types (NOQT) to be 



FQH state: 




$3. 7 






$6.7 




-,z s 


$9. 7 


$10. 7 

— ;^io 


NOQT: 


3 


4 


5 


6 


7 


8 


9 


10 


11 



For the parafermion states $ _^_ ;Z „ = $f ; z„ IIi<j( z i — z j) 2 (™ = 2 + 2n), we find 



FQH state: 


$2. 7 


$3.7 


$4.7 


$ « -z. 


$6.7 


$7.7 


$8.7 


$9.7 


$10.7 


NOQT: 


9 


16 


25 


36 


49 


64 


81 


100 


121 



For the generalized parafermion states $„_. z (no, we find 



FQH state: 


$5. z <2) 

£J J 


$ 5 . z (2) 


$7.7(2) 
R > Z 7 


$ 7 .7(2) 
55 >^7 


$ 7 , 7 (3) 
TS 1^7 


$ 7 . z (3) 
35 '^7 


$ 8 . z (3) 


$9. 7 <2) 
R 


$10.7(3) 

1R>^10 


NOQT: 


24 


54 


32 


88 


72 


128 


81 


40 


99 



where fc and n are coprimc. For the composite parafermion states $ ni „(k 2 )$„ 2 „(k 2 ) obtained as products of two 

— ^"i — ;^«2 

parafermion wave functions , we find 



FQH state: 




*§iZ,*i;Z. 


$ |;Z2 $ |;Z 5 




NOQT: 


30 


70 


63 


117 



where n\ and n 2 are coprime. The filling fractions of the above composite states are v = 



The above results suggest a pattern. For a (general- 
ized) parafermion state $ J£L . z (*o, we can express its fill- 
ing fraction as v = n/m = p/q where p and q are co- 
prime. Then the number of quasiparticle types is given 
by NOQT= qD{n) where D(2) = 3, D(3) = 2, D{A) = 5, 
D(5) = 3, D(6) = 7, D(7) = 4, D(8) = 9, D(9) = 5, 
and D(1Q) = 11; or D(n) = n + 1 for n = even and 
D(n) — ^y^- for n — odd. Similarly, For a composite 
parafermion state $ a (^i^^ (t 2 ), we can express its 

m 1 ' "1 m 2 ' "2 

filling fraction as ^ = ^i 23 = p/q where » and q 

min 2 +m 2 ni ' ^ ^ 

are coprime. Then the number of quasiparticle types is 
given by NOQT= g£>(ni)D(n 2 ). 

The corresponding CFT of the above (generalized and 
composite) parafermion states are known. The numbers 
of the quasiparticle types can also be calculated from the 



CFT. 4 For the generalized parafermion state $n_. z (*) 
the numbers for the quasiparticle types is given by 



34 



NOQT = 1 n ( n+1 ) = m n{n+l) = m(n + 1) 
V v 2 n 2 2 

For the composite parafermion state Ili $_2±..z ( ' c i ) * nc 
numbers for the quasiparticle types is given by 34 

NOQT= l -\\ n ^ + l) . 
v - L - L 2 

i 

Here we require that fcj is not a factor of and m, ri2, TI3, 
■ ■ ■ have no common factor. The CFT approach gives rise 
to exactly the same numbers for the quasiparticle types. 



For generalized parafermion states $ n.. z (*) where n and k have a common factor, we have 



FQH state: 


$4.7(2) 


$6.7(2) 


$6. 7 (3) 
R > Z 6 


$8.7(2) 


$8.^(4) 
R.^R 


$ 9 .7(3) 


NOQT: 


10 


20 


21 


35 


36 


56 



For more general composite parafermion states $ ni z (k 2 )^ n 2 .^fe) where n\ and n 2 have a common factor, we have 

m 1 > "1 m 2 ' "2 



FQH state: 


$2.7 $2.7 

2 ;^2 2 ;-Z 2 


$2.7 $4.7 

2 ,z 2 2 ,z 4 


$3.7 $3.7 
2 ,z 3 2 ,z 3 


$4.7 $4.7 

2 ,^4 3'^ 4 


$5. 7 $5.7- 




$5. 7(2) $5. 7(2) 
8^5 5^5 


NOQT: 


10 


42 


20 


35 


56 


352 


224 



For the v = 1/2 states $^. ; c n , we have 



FQH state: 


$3.^ 

6>C3 




$ 5 - 

To.L-s 


$ 6 - 


$ 7 -r, 


$ 8 .~ 


$ 9 ~ 


NOQT: 


56 


170 


352 


910 


1612 


3546 


6266 
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Note that ^^_. C r an d $s. z $ s < 2 ) have the same pattern of zeros and may be the same state. For the v — 1 states 
$£;C„/z 2 ! we nave 



FQH state: 


$ !;C 4 /Z2 




$ l;Cs/Z2 




NOQT: 


35 


138 


171 


338 



Note that <E> 8 _r 2 ), $4.^ , 7 and $4. 7 <I>4. 7 have the same pattern of zeros and may be the same state. 

%;Z\ " 4,04/^2 2^4 2< z 4 



For those more general and new FQH states, the cor- 
responding CFT are not identified. Even the stability of 
those FQH states is unclear. If some of those states con- 
tain gapless excitations, then the number of quasiparticle 
types will make no sense for those gapless states. In the 
following, we will study a few simple examples in more 
details. 



B. v = 1/2 Laughlin state 

For the v = 1/2 Laughlin state, n — 1 and its pattern 
of zeros is characterized by 



$1/2 : (m;S 2 , ■ 
(-■h sc ■ 
n ■ 



■ , S n ) — (2; ), 

•,C) = 
■n m -\ = 10 



(f;0), 



Solving (22), we find that there are two types of quasi- 
particles. Their canonical occupation distributions and 
other quantum numbers are given by 









10 
01 




1/2 



1/4 



Thus Ti, T 2 , and T are given by 
Ti = 






Since l max = 0, from (92), we find that the conjugate of 
10 is TO = 10 and the conjugate of 01 is UT = 01. Thus 
C= 1. 

From TT 2 = c m l m T 2 T x T (see eqn. (49)) we find that 
ti = i. From STi = T 2 S, we find that 5i = #00 and 
S01 = -Sn. From S 2 = C = 1, we find that S n = -S 00 , 



and 5*1 



00 



± 



V2' 



Thus 





T 



The above implies (ST) 3 — e 17r / 4 . Those modular trans- 
formations, S and T, agree with those calculated using 
the Chern-Simons effective theory 20 and the explicit FQH 



wave functions. 36 From Tn = i, we find that the quasi- 
particle (rij-i) = (01) has a scaling dimension ft = 1/4 
and a semion statistics. 
Let us introduce 

rp _ i2n/2irp c _ c 

We find that T m and S m satisfy 



1. 



( S m T m ) 



1. 



Thus S m and T m generate a linear representation of the 
modular group. 



C. Z2 parafermion state 

The v = 1 bosonic Pfaffian state 35 is a Z 2 parafermion 
state with n = 2. Its pattern of zeros is described by 

$| ;Z2 : (m;5 2 ,--- , S n ) - (2;0), 

(^;ftf,...,ft- ) = ( 2.1 ) 0), 

(n ,-- • ,n m _i) = (2,0) 

Solving (22) , we find that there are three types of quasi- 
particles: 



n T,l 


Q-f 


5° s p 


20 








02 


1 


1/2 


11 


1/2 


1/4 





Thus in the I7) basis, T\, T 2 and T are given by 



To = 



SmceJ_ max = 0, from (92), we find that 20 = 20, 02 = 02, 
and TT= 11. Thus C = 1. 

From TT 2 = e i ™/ m T 2 T 1 T (see eqn. (49)) we find that 
t\ = —1, but t 2 is undetermined. Using ST\ = T 2 S, 
ST 2 = TxS, S T = S, and S 2 = 1, we obtain the following 
possible solutions: 



S = 



1 



' 1 1 V2x 

1 1 -V2a; 

v V2/a; -V2/x 
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Those S's satisfy (ST) 3 = t 2 . 

Using the above S, we can calculate the fusion coeffi- 
cients <5 7l7273 and S7T72 from eqn. (71). We find that 
6220 = S221 = 1/x 2 . The condition (5 77 *o = 1 fixes 
x = ±1. Thus we have the following two possible so- 
lutions: 



Si = - 



and 



S 2 = 



We note that the above 5"s are already symmetric. Thus 
those S"s can regarded as S TC . 

Now let us calculate A^ and 7 = 0,1,2 using 

eqn. (70). For the first solution S\ we find that 



(4 7 °) 



where 7 labels rows and 7' labels columns. We have 




A (0) 



A (2) = 



We see that 




where I7), 7 = 0, 1, 2 are the degenerate ground states. 
For the solution S 2 we have 



(4 y) ) 



(1 1 


V2) 




1 1 


-V2 


and 


\V2 -y/2 


J 




( 1 1 


-V2] 




1 1 






\-\/2 V2 







( A (o) A (o) aWaw a^aw\ ( aw aw aw 
AW AM AW AW AW AW = AW AW AW 
A (2) A (0) A (2) A (1) A (2) A (2)J ^(2) A (2) A (0) +A (l) 

We recover the fusion algebra of the Z2 parafermion the- 
ory. We also find that 

(l o\ 

B (o) = 



BW 




bW = 



/o 1 o\ 

1 

V° V 



We see that B^- 1 ' also encode the fusion algebra 

(\0) |1> 
|1> |0) 



/b(0)| ) B(0)|l) fi(0)|2) ' 
SW|0) Bf 1 )!!) Bf 1 )^) 
V B(2)| ) B( 2 )|l) BW\2) y 



|2> 
|2> 

VI2) 12) |0> + |1>, 




A (0) 



A (2) = 




/l 
1 

vQ -1, 



BW) remain the same. 

The first row of the 5-matrix are called quantum di- 
mensions of the quasi-particles. For a unitary topological 
quantum field theory, all quantum dimensions must be 
positive real numbers, and moreover > 1. Therefore, the 
second solution does not give rise to a unitary topologi- 
cal field theory. Based on this reason, we exclude the S2 
solution. 

The v = 1 bosonic Z2 parafermion state &i-z 2 nas 
three degenerate ground states on a torus. In the thin 
torus limit, the three ground states are described by the 
occupation distributions 20202020 • • • , 02020202 • • • , and 
11111111 • • •. 

The v = 1/2 fcrmionic Z 2 parafermion state 
Y[( z i ~ z j) nas srx degenerate ground states on a 
torus. In the thin torus limit, the six ground states are 
also described by the occupation distributions. Those 
occupation distributions can be obtained from that of 

I$2. Z2 state given above. We note that multiplying the 
factor n( z »~ z j) increases the space between every neigh- 
boring particles in a distribution by one. For example, 
it changes 11 to 101, 101 to 1001, 2 to 11, 3 to 111, 
etc . It changes the bosonic distribution 202020 • • • to 
a fermionc distribution 110011001100 ••• , and changes 
the distribution 111111 ••• to 101010101010 ••• . In- 
cluding the translated distributions of 110011001100 • • • 
and 101010101010- ■-, we find that the fermionic Z 2 
parafermion state 3>| ; z 2 Y\(zi — zj) has six degenerate 
ground states described by the distributions 

110011001100 
011001100110 
001100110011 
100110011001 
101010101010 
010101010101 

Note that a unit cell contains to = 4 (or 2) orbitals. 
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D. Z3 parafermion state 

The v = 3/2 bosonic Z 3 parafermion state has a pat- 
tern of zeros described by 

$s. Z3 : {m;S 2 ,--- ,S n ) = (2; 0,0), 



(— \h 

n 



l ' 

(n ,- • • ,n m _! 



2 2 2 
,^ = (3; 3, g,0), 



(3,0) 



Solving (22), we find that there are four types of quasi- 
particles: 





O7 




30 








03 


3/2 


3/4 


12 


1 


1/2 


21 


1/2 


1/4 



Thus in the I7) basis, Ti, T 2 and T are given by 



Ti = 



/l 





\ 





-1 











1 


Vo 





-1) 


/1 





\ 





h 











t 


\° 





htj 



To = 



/o 1 o\ 

10 
1 
\0 1 0/ 



Since l max = O^from (92), we find that 30 = 30, 03 = 03, 
21 = 21, and 12 = 12. Thus C = 1. 

From TT 2 = c m l m T 2 T{T (sec cqn. (49)) wc find that 
£1 = — i and t 2 = t. Thus 



T - 



(l \ 

-i 

i 

\o — i*y 



We can rewrite Ti, T2, and T in direct product form 




Using STi = XaS, ST 2 = T X S, and 5 2 = 1, wc find that 
S must have the following form 



cos(0) sin(0) 
sin(0) -cos(0) 



Note that the above S is already symmetric. Thus we 
can regard the above S as the S TC . 




The direct product form of Ti, T 2 , T, and S sug- 
gests that the U(l) charge part and the non-Abelian part 
separate. 28 So let us concentrate on the non-Abelian part 



T = 




cos(0) sin(0) 
sin(0) — cos(0) 



From S, we can calculate and [see eqn. (70) 

and eqn. (69)]. All the matrix elements must be 

non-negative integers and (ST) 3 oc 1. 

One way to satisfy those conditions is to let 



tan(0) 



In this case 



T = 




1 + V5 



S = 




(74) 



and 



A m = hp 

10 l-tp. 




A<n = 1, 
P = 1, 



Those are the only valid solutions, which are realized 
by the Fibonacci anyons. The solutions with d\ < are 
excluded as reasoned at the end of the last section, and 
one of the them is the non-unitary Yang-Lee model. 

Taking 71 = 72 = 1 in cqn. (72), we find that 



-1 



or 



cos(2ttIii 



-<P/2, 



which gives us 



hi = ±2/5 mod 1. 



Since ho = 0, we see that e l27rh ~< are the eigenvalues of the 
T operator (see eqn. (74)). Let us denote the eigenvalues 



of the T operator as e 



Wc find, for the two choices 



of h{ = ±2/5 (see eqn. (74)), 





h i 




30 








03 


3/4 


3/4 


12 


2/5 


3/5 


21 


3/20 


7/20 



According to CFT (see appendix 6), the scaling dimen- 
sion of the quasiparticlc operators in the Z3 parafermion 

Q 2 

FQH state is h~ t — + h^. For the quasiparticles 30 



19 



and 03, h 1 ^ = since those are Abelian quasiparticles. 
For the quasiparticles 12 and 21, = 1/15. Thus 
h n = 0, 3/4, 2/5, and 3/20 for the quasiparticles 03, 
30,12, and 21 respectively, which exactly agree with 

for the case of hj = 2/5. This example demonstrates a 
way to calculate 27 quasiparticle scaling dimensions from 
the pattern of zeros. 

The v = 3/2 bosonic Z3 parafermion state ^| ; z 3 nas 
four degenerate ground states on a torus, described by 
the occupation distributions 

3030303030 • • • 
0303030303 • • • 
2121212121 • 
1212121212 • 

The v = 3/5 fermionc Z 3 parafermion state $|. Z3 n( 2 i — 
Zj) has 10 degenerate ground states on a torus, described 
by the occupation distributions 

1110011100111001110011100 • 
0111001110011100111001110 • 
0011100111001110011100111 • 
1001110011100111001110011 • 
1100111001110011100111001 • 
1101011010110101101011010 • 
0110101101011010110101101 • 
1011010110101101011010110 • 
0101101011010110101101011 • 
1010110101101011010110101 • 
Note that a unit cell contains m = 5 orbitals. 
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APPENDIX 
1. The calculation of orbital spin 

F L ,M) = - J)c- a2[2 S l%m - L 

1=0 

= Y,( km + L -J)^ a2(km+L)2 - (75) 



Let 



k=0 



We find that, for a canonical occupation distribution 
(25), we can rewrite (31) as 



m— 1 



6l = nY,{n T ,i-ni)F l {a). 



(76) 



To evaluate Fl,j(cx), we will use the Euler-Maclaurin for- 
mula: 

1 k=i J o 

= E(^(/ w (^)-/ w (o))+fi, 



where Bk are Bernoulli numbers (B , B\, B 2 , B 3 , ■ ■ ■ ) 
.1 I — — 

2' 6> ' 30' 



VI. CONCLUSIONS 

Through string- net wave functions 5 ' 6 , one can show 
that non-chiral topological orders can be naturally de- 
scribed and classified by tensor category theory. This 
raises a question: how to describe and classify the chiral 
topological order in FQH states? The results in Ref. 11 
suggest that the pattern of zeros may provide a way to 
characterize and classify chiral topological orders in FQH 
states. In this paper, we see that many topological prop- 
erties of chiral topological orders can be calculated from 
the data {S a } that describe the pattern of zeros. In par- 
ticular, through the algebra of tunneling operators, we 
see a close connection to tensor category theory. The pat- 
tern of zeros provides a link from electron wave functions 
(the symmetric polynomials) to tensor category theory 
and the corresponding chiral topological orders. 
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\R\ < 



(^) 2p Jo 



\f^ +1 Hx)\dx. 



For our case f(x) = (L + mx — jj c - a2 ( mx + L ) 2 _ if wc 
choose a ~ 1/J and p = 1, in large J limit, we find 
R ~ 1/J in large J limit. Therefore 



F L A a ) 

= I (L + mx- J)e- a2{mx+L ^ Ax 
Jo 

+ \(L-J)-^+0(l/J) 



OO J I X T 



ma 



\(L-J)-^ + 0(l/J) 



f 00 L + — — J 

\ -f (l-2aLx- a 2 L 2 + 2a 2 L 2 x 2 )c- x dx 

Jo 



ma 



+ \{L-J)-^ + 0(l/J) 
olJ\Fk -1 J r L 2 L - J m 

v I T I 

2a 2 m m 2m 2 12 



0(1/ J). 
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a unit cell 



12 /, 



FIG. 11: The graphic representation of an occupation dis- 
tribution (n ,--- , nio) = (2,0,1,0,2,0,0,0) for the <E> 5 „(a) 



state with n — 5 and m = 8. L 



; 4 for such a distribution. 



Since X^lio ( n ' — n 7;0 = for a canonical occupation 
distribution, the terms that do not depend on L do not 
contribute to the Berry's phase. Thus we have 



m j 1 / 2 

0l = nJ2(n r , l -n l )(-l + ---). (77) 

1=0 



Compare with (26) and (28), we find that the -^l term 
exactly reproduces the quasiparticle charge. The terms 
that do not depend on J give us the orbital spin (32). 

Since ni is a periodic function of I with a period m, 
we may also view n; as a periodic function with a period 
km. If we view m as a periodic function with a period 
/cm, the orbital spin S° sp will be given by 



km — 1 , .2 

1=0 



7 



2 2fcm 

jm + I (Jm + I) 2 



2km 



k-l m-l 

3=0 1=0 
k— 1 m— 1 

]=0 1=0 

./kl k(k-l)l I 2 \ 
= 2^( n 7;i-^)(- ^ — ) 



2km 



1=0 
m— 1 

E 

(=0 



(n 7;i - n^f 



2k 

2 2m) 1 



2m) 



which is identical to the previous result (32). Therefore, 
the formal calculation of S° sp using e~ a 1 regulator (see 
(31)) produces a sensible result. 



2. Orbital spin of Abelian quasiparticles 

If a quasiparticle is described by an occupation dis- 
tribution n T j that can be obtained by shifting the oc- 
cupation distribution ni for the ground state, then the 
orbital spin of such a quasiparticle can be calculated re- 
liably without using the formal unreliable approach de- 
scribed above. Such kind of quasiparticles can be created 
by threading magnetic flux lines through the FQH liquid 
and correspond to Abelian quasiparticles. 

The occupation distribution n; for the ground state has 
some properties that will be important for the following 




1 2 



FIG. 12: If JV^, — U — lmax+N c m, the symmetric polynomial 
$ of N = N c n variables can be put on the sphere without any 
defects (or quasiparticles). 



I I 

i i 3 n 3 i i 3 n 3 



1 2 



2.1 



FIG. 13: If = 2 J = l m ax + (N c - l)m + n^, the above 
occupation distribution n 7; ; (obtained by shifting the ground 
state distribution in Fig. 12) describes a single quasiparticle 
on the sphere located at the south pole (1 = 0). There is no 
quasiparticle at the north pole (I = 2J). 



discussion. In addition to the periodic property ri; +r , 
ni, ni also have a symmetric property 



ni = ni 



0<l< L 



(78) 



according to numerical experiments, where l ma x = S n — 
S n -i is the largest I in the first unit cell < I < m such 
that ni max > (see Fig. 11). ((78) implies that no > 0). 

On a sphere with = 2 J flux quanta, there are A^+l 
orbitals labeled by I = 0, 1,- • • , 2 J (see Fig. 2). Those 
orbitals form a angular momentum J representation of 
the 0(3) rotation. The J z quantum numbers of those 
orbitals are given by m z = I — J. If the flux through 
the sphere is such that = 2J = l max + (N c — l)m 
for an integer -ZV C , then the occupation ni can fit into 
the 2 J + 1 orbital in such a way that there is no non- 
trivial quasiparticle at the north and south poles (see Fig. 
(12)). Such a state has N = N c n particles. Such an N c n- 
particle state can fill the sphere without any defect and 
form a J tot = state. Note that l max = l n — S n — S n -x 
and S n _i = (see Ref. 11), thus l max = ^S n . So an 

iV c n-particle state can fill the sphere without any defect 
if Nj, — 2J = 2-S n + (N c — l)m which is exactly the 
condition obtained in Ref. 11. 

Let us create a quasiparticle by threading flux lines 
through the south pole. The total flux quanta becomes 
— 2 J — ^S n + (N c — l)m + and the occupation 
distribution n 7; ; for the created quasiparticle is obtained 
from m by shifting the distribution by (see Fig. 13). 
The occupation distribution is identical to that of the 
ground state distribution in Fig. 12 near the north pole 
(/ = 2 J). Thus the distribution describe a state that 
has no quasiparticle near the north pole. However, the 
occupation distribution is different from the ground state 
distribution near the south pole (I = 0). Therefore, the 
distribution describes a state with a quasiparticle near 
the south pole. 

The total J z of the above quasiparticle state is given 
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by 



2J 



2J 



1=0 1=0 
_ n^n / TO - Irnax _ n <t>\ 

~ m V + 2 2) 



JQ 1 + ?L-^Q y -HLQ2 



where 



(79) 



is the quasiparticle charge. If we move such quasiparti- 
cle along a loop that spans a solid angle SI, the induced 
Berry's phase (P B will be 

®B _ Ttot If^i i m ~ ^max n _ m n 2 

fi ~ z ~ ^ 7 2 ^ 7 2n^ r 

Compare with (30), we find the orbital spin of the quasi- 
particle to be 



~ 2 Ql 2n^ 



7 



(80) 



Let us compare (80) with (32). For the Abelian quasi- 
particle, its occupation distribution n 7; ; has a form 

In this case, (32) becomes 



S os P = J2 ni (^ l + n ^ _ ( l + n <p) 2 _L + Jl_\ 
i=o 2 



2m 2 2m) 

m— 1 



2m ^ m 

1=0 

Using (78) , we find that the above expression agrees with 
(80). This confirms the validity of (32) and (33) for the 
case of Abelian quasiparticles. On the other hand, the 
validity of (32) or (33) for the case of non- Abelian quasi- 
particle is yet to be confirmed by a more rigorous calcu- 
lation. 



3. Translations of orbitals on torus 

To show T 2 4> {1) = (f> {l+1 \ we note that 
T 2 ^ l \x\X 2 ) 

1 A i7 rAVr(X 2 + 1 i-) 2 



c i2 ^ 1 / ( ' ) (^ 1 +r(A 2 + ^))c i 



Since 



i w- (N^k+l) 2 + i 2ir( JV k+l) (z+ J- ) 



i^-(N 4 ,k+l+l) 2 + i2TT(N 4 ,k+l)z 



k 



we have 

T 2 ^ l \x\X 2 ) 

= e i27rXl - i ^- i2w 7( i + 1 )(A: 1 +TX 2 )e i ^ r(x2+ ^ )2 
= <f>V +1 \X\X 2 ). 



Non-Abelian Berry's phase and modular 
transformation 



The wave functions <j>{„^. i+1 j form a basis of the de- 
generate ground states. As we change the mass matrix 
or r, we obtain a family of basis parameterized by r. 
The family of basis can give rise to non-Abelian Berry's 
phase 37 which contain a lot information on topological 
order in the FQH state. In the following, we will discuss 
such a non-Abelian Berry's phase in a general setting. 
We will use 7 to label the degenerate ground states. 

To find the non-Abelian Berry's phase, let us first de- 
fine parallel transportation of a basis. Consider a path 
t(s) that deform the inverse-mass-matrix g{r\) to g(r 2 ): 
T\ = t(0) and t 2 = r(l). Assume that for each inverse- 
mass-matrix <7[t(s)], the many-electron Hamiltonian on 
torus (X\X 2 ) ~ (A 1 + 1,A 2 ) ~ (A 1 , A 2 + 1) has N q - 
fold degenerate ground states I7; s), 7 = 0, 1, • • • , N q — 1, 
and a finite energy gap for excitations above the ground 
states. We can always choose a basis I7; s) for the ground 
states such that the basis for different s satisfy 



(7';s| — |7;s) = 0. 

Such a choice of basis I7; s) defines a parallel transporta- 
tion from the bases for inverse-mass- matrix g(ri) to that 
for inverse-mass-matrix g(r 2 ) along the path r(s). 

In general, the parallel transportation is path depen- 
dent. If we choose another path r'(s) that connect n 
and t 2 , the parallel transportation of the same basis for 
inverse-mass-matrix g(r\), I7; s = 0) = \y;s = 0)', may 
result in a different basis for inverse-mass-matrix g(r 2 ), 
|7;s = 1) 7^ \ f,s = 1)'. The different basis are related 
by an invertible transformation. Such a path depen- 
dent invertible transformation is the non-Abelian Berry's 
phase. 37 

However, for the degenerate ground states of a topolog- 
ically ordered state (including a FQH state), the parallel 
transportation has a special property that, up to a total 
phase, it is path independent (in the thermal dynamical 
limit). The parallel transportations along different paths 
connecting t\ and t 2 will change a basis for inverse-mass- 
matrix g{r\) to the same basis for inverse-mass-matrix 
t(t 2 ) up to an overall phase: I7; s = 1) = e 1( ?|7;s = 
1)'. In particular, if we deform an inverse-mass-matrix 
through a loop into itself (ie r(0) = t(1)), the basis 
I7; 0) will parallel transport into I7; 1) = e iip \j: 0). Thus, 
non-Abelian Berry's phases for the degenerate states of 
a topologically ordered state are only path-dependent 
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Abelian phases e 1(p which do not contain much infor- 
mation of topological order. 

However, there is a class of special pathes which give 
rise to non-trivial non- Abelian Berry's phases. First we 
note that the torus (X\A 2 ) ~ (X 1 + 1,X 2 ) ~ (A\A 2 + 
1) can be parameterized by another set of coordinates 





X ' x' x'l o x' 

(a) (b) (c) (d) 



FIG. 14: (a): A system described by an inverse-mass- matrix 
g(r). (b): A system described by an inverse-mass-matrix 
g(r') where r and r' are related by a modular transforma- 
tion Mr- (c): The system in (b) is the same as the system in 
(a) if we make a change in coordinates, (d): The system in 
(c) is drawn differently. 



where a,b,c,d e Z, ad — be = 1. The above can be 
rewritten in vector form 

X' = M~ X X, X = MX', M e SL{2, Z). (82) 

(X n ,X' 2 ) has the same periodicity condition 
(X'\A' 2 ) ~ (X n + 1,X' 2 ) ~ (X n ,X' 2 + 1) as 
that for (A 1 , A 2 ). We note that 



-b a I \d x ,2 \d X '2 \b d \d X 2 



The inverse-mass-matrix in the (A" 1 , A" 2 ) coordinate, 
g(r), is changed to 





in the (A' 1 , A' 2 ) coordinate. From eqn. (36), we find 
that 



= 9{r'), 



with 



b + dr 
a + ct 



(83) 



(84) 



The above transformation t — > r' is the modular trans- 
formation. We see that if t and r' are related by the 
modular transformation, then two inverse-mass-matrices 
g(r) and g(r') will actually describe the same system 
(upto a coordinate transformation). 
Let us assume that the path t(s) connects two r's 

(:;> 



related by a modular transformation M 



r(l) 



^ct^) ■ We WU1 denOt0 T ( ) = T 2111(1 T (!) = T '' 

The parallel transportation of the basis I7; r) for inverse- 
mass-matrix <?[t(0)] gives us a basis |7;r') for inverse- 
mass-matrix gr[r(l)]. Since r = r(0) and t' = t(1) arc 
related by a modular transformation, <?[t(0)] and ff[r(l)] 
actually describe the same system. The two basis I7; r) 
and 1 7; t') are actually two basis of same space of the 



degenerate ground states. Thus there is an invertible 
matrix that relate the two basis 



\ r ,r')=U(M)\r,r) 
C/ 7 y(M) = (7;r|7 / ;r / ) = (7; t|*7(M)| 7 '; r). 



(85) 



Such an invertible matrix is the non-Abelian Berry's 
phase for the path r(s). Except for its overall phase 
(which is path dependent), the invertible matrix U is a 
function of the modular transformation M. In fact, the 
'invertible matrix U form a projective representation of 
the modular transformation. The projective represen- 
tation of the modular transformation contains a lot of 
information of the underlying topological order. 

Let us examine (7;t|7';t) in eqn. (85) more care- 
fully. Let $ 7 [{AQ}|t] be ground state wave functions 
for inverse-mass-matrix g [r], and $ 7 [{Xi}|T'] be ground 
state wave functions for inverse-mass-matrix g[r']. Here 
Xi = (A*,A 2 ) are the coordinates of the i th electron. 
Since r and r' are related by a modular transformation, 
$ 7 [{Xi}|r] and $ 7 [{Xi}|r'] are ground state wave func- 
tion of the same system. However, we cannot directly 
compare $ 7 [{Xi}|r] and $ 7 [{Xi}|T'] and calculate the 
inner product between the two wave functions as 

($ 7 [{XJ|t(0)], <M{XJ|t']). 

The wave function $ 7 [{Xi}|r'] for inverse-mass-matrix 
g[r'] can be viewed as the ground state wave func- 
tion for inverse-mass- matrix g[r] only after a coordi- 
nate transformation (see Fig. 14). Let us rename 
X to X' and rewrite $ 7 [{XJ|t'] as $ 7 ({A^}|r'). 
Since the coordinate transformation (82) change r 
to t', we see that we should really compare 
$ 7 ({X'J|t') = $ 7 ({M" 1 XJ|r') with $ 7 ({X 4 }|r). But 
even <f> 7 ({Xi}|r) and <i> 7 ({M _1 Xj}|r') cannot be di- 
rectly compared. This is because the coordinate trans- 
formation (81) changes the gauge potential (37) to an- 
other gauge equivalent form. We need to perform a U(l) 
gauge transformation Ug(M) to transform the changed 
gauge potential back to its original form eqn. (37). So 
only $ 7 [{XJ|t] and U G ^-y[{M- 1 X i }\T'] can be directly 
compared. Therefore, wc have 

C/ 77 ,(M) (86) 

= ($ 7 ({Xj|r), £/ G (M)<iy ({M^XOIr')) 
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which is eqn. (85) in wave function form. Note that r = 

M12+M221 
M11+M21T ' 

Let us calculate the gauge transformation Ug(M). We 
note that <& 7 ({X-}|t') is the ground state of 



k QY ,Z 

k " »,j=l,2 fe 



ax? J 



where fc = l,--- ,7V labels the different electrons. In 
terms of X 1 (sec cqn. (81)), H' has a form 



4^ 2 . f-' 'Mi 



fe " i,j=i,2 N ^ 



- \A 



where 




Since (A' 1; ^) = (-27r/V X' 2 , 0), Wc find that 

(A 1 ,A 2 ) = {-2irN^X' 2 d, 2nN^X' 2 b) 
= {^irN^-cX 1 + aX 2 )d, 2irN (t> (-cX 1 + aX 2 )b) 

U G (M) will change H' to H: 



U G (M)H'U G (M) = H 



d 



with (^1, A 2 ) = (-27r7V A 2 ,O). Wc find that 

C/ G (M) = ]J WG (M;X fe ), 

k 

u G {M;X) = j**N+[bcX*x*-%W-$(x>n (87) 

Eqn. (86) can also be rewritten as a transformation on 
the wave function $ 7 ({Xi}|r): 

U(M)^({Xi}\ T ) = U G {M)%{{X' i }\r') 

= t^ G (M)* 7 ({M- 1 X i }|T y ) 

= $ r ({X i }|r)f/y 7 (M). (88) 

where r' = ^ 12 ~l~ w 22T . We see that the action of the op- 
erator U(M) on a wave function $ 7 ({Xi}|r) is to replace 
X t by X\ = M~ 1 X l , replace r by r' = $Lt%toz , an d 
then multiply a phase factor U G (M) given in eqn. (87). 
Thus eqn. (88) defines a way how modular transforma- 
tions act on functions. We find that 



U(M'){U{M)^{{Xi}\T) 
= U(M')(u G (M)^({M- 1 X i }\r') 

= U G (M' M)$ 7 ({M" 1 M'" 1 Xi\\T lr 
= C/ g (M'M)$ 7 ({(M'M)- 1 X 4 }|t") 
= U(M'M)^({Xi}\T). 



nere r - Mii+MaiT and t - Mll +M 21 f wltn T _ 

U(M')U(M) = U(M'M). 

So U (M) form a faithful representation of modular trans- 
formations SL(2, Z). 

To summarize, there are two kinds of deformation 
loops t(s). If t(0) = t(1), the deformation loop is con- 
tractible [ie we can deform the loop to a point, or in 
other words we can continuously deform the function 
t(s) to a constant function r(s) = t(0) = r(l)]. For 
a contractible loop, the associated non-Abclian Berry's 
phase is actually a U(l) phase U 11 i = e ltp <5 77 '. where 
if is path dependent. If r(0) and t(1) are related by 
a modular transformation, the deformation loop is non- 
contractiblc. Then the associated non-Abelian Berry's 
phase is non-trivial. If two non-contractible loops can 
be deformed into each other continuously, then the two 
loops only differ by a contractible loop. The associated 
non-Abelian Berry's phases will only differ by an over- 
all U(l) phase. Thus, upto an overall U(l) phase, the 
non-Abelian Berry's phases C/ 77 ' of a topologically or- 
dered state are determined by the modular transforma- 
tion t — ► t 1 = ~~T~j • We also show that we can use 
the parallel transportation to defined a system of basis 
$ 7 ({X i }|r) for all invcrse-mass-matrices labeled by r. 
By considering the relation of those basis for two t's re- 
lated by an modular transformation, we can even obtain 
a faithful representation of the modular transformation 
SL(2,Z). 



5. Algebra of modular transformations and 
translations 

The translation Td and modular transformation U(M) 
all act within the space of degenerate ground states. 
There is an algebraic relation between those operators. 
From (88), we see that 



d}\r) 



[/(M)T d $ 7 ({XJ|r) = *7(M)$ 7 ({X, 
= U G $ 7 ({M- 1 X i + d}\T') 
= U G <S> 1 {{M- 1 {X l + Md)}\ T '). 
Therefore 

e ie U(M)T d = T Md U(M). (89) 

Let us determine the possible phase factor e'^ for some 
special cases. Consider the modular transformation r — > 
t' = t + 1 generated by 



M T = 



1 1 



M; = 



1 -1 



11 1 V° 1 / 

We first calculate 

^(M T )Ti$ 7 ({X<}|r) = U{M T )<f> 7 {{Xi + di}\r) 
= e- i7rAr * s ^ ( ^ 2)2 $ 7 ({M T 1 X 4 + d 1 }|r'). 
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where di = (-^-,0). We note that 



Thus we next calculate 

T!;7(m t )$ 7 ({xj|t) 

^T 1 c-^ n ^ {x ^ 2 ^{{M t 1 X 1 }\t') 

= e- i " Ar *^ ( ^ 2)2 $ 7 ({M T 1 (X l + d 1 )}\r') 

= c -^ Ar *^( x ') 2 $ 7 ({M T 1 X, + d 1 }|r'). 

Therefore 

U{M T )T 1 = TiC/(M T ). 

To obtain the algebra between U{Mt) and T 2 , we first 
calculate 

U(M T )T 2 ^({Xi}\T) 
^U{M T )e^^ x ^ 1 {{X t + d 2 }\r) 
= e-^ N ^ x ^\ i ^^ x '- x ^%{{M T 1 X l + d 2 }\ T '), 

where d 2 = (0, ^-). We note that 



MTd2 = o 1 



Thus we next calculate 

T 2 7W(M T )$ 7 ({X,}|t) 
= T 2 T 1 c-' 1 * N f^ x ^ 2 <5> 1 {{M T 1 X l }\T') 
= T 2C -^ N ^ x2s > 2 ^ 1 {{M T \X l + d 1 )}\T l ) 

$ 7 ({Mf 1 (A: j + di + d 2 )}|7 y ) 



We see that 



U{M T )T 2 = c i7r ^T 2 7W(M T ). 

Next we consider the modular transformation r 
t' = — 1/t generated by 



M s = 



-1 
.1 



Ms 1 



o r 
-i o, 



From 



we find 

[/(M 5 )[/(M 5 )$ 7 ({XJ|t) 

= U{M s )e- i2 * N ^ x ' x2 %{{Ms 1 X i }\T') 

= e -i2^ £, x, 1 *? e -i2^, E, ^?(-^* 1 )$ 7 ({_x i }|T) 

= $ 7 ({-Xj|r). 

We see that U(M S )U(M S ) = U(-l) generates to trans- 
formation X — > —X. We can show that 

L^-lJTiCT^-l) = Tf 1 , [/(-^TsCZ-^-l) = Tjj- 1 . 

(90) 

Since the wave function of an orbital satisfies 

<t>«\-x\-x 2 ) = ^ N *- l \x\x 2 ) = ^- l \x\x 2 ), 

we find 

U(M s )U{M s )%{{Xi}\T) = ^(-l)$ 7 ({Xi}|r) 

= * 7 .({X i }|r), (91) 
where 7* corresponds to the occupation distribution 

Since n 7i ; and n 7 *^ are periodic with a period of m, the 
above can be rewritten as 

n Y;i = n T,(l max -l)%m- (92) 

7* also corresponds to the anti-quasiparticle of 7. Thus 
we find that 



[U(Ms)} 2 = C, 



C$ 7 = $„„. 



[/(M s )$ 7 ({X,}|r) - c- i2 ^^ x . lx ?$ 7 ({M s 1 XJ|T'), 



where C is the quasiparticle conjugation operator. 
Clearly C 2 = 1. 
Let us first calculate 

C/(M s )T 1 $ 7 ({XJ|r) 

= C/(A/f s )$ 7 ({M s 1 X t + d 1 }|r) 

= e- i2wJV *^ x . lx . 2 $ 7 ({M^ 1 Xi + d 1 }|T / ). 
Since M$d\ = d 2 , we next consider 
T 2 C/(M s )$ 7 ({XJ|r) 
= T2 e- i2 ^E^, 1 ^ 2 $ 7 ({ M -i Xi }| T / ) 

^({M^^Xi + da)}!/) 

= e - i27rJV *^ x « lx ?$ 7 ({M5 1 (Xi + d 2 )}|T / ). 
We see that 

U(M S )T! = T 2 U{M S ), 

U(M S )T 2 = T^U(Ms), 

where we have used (90). Let us introduce 

T=e i9 U(M T ), S = U(M S ), 

where the value of 8 will be chosen to make T 00 = 1 . We 
find that 

TT\ = T\T, TT 2 = e i7r -T 2 TiT, 

STx = T 2 S, ST 2 = CTiC~ x S. (93) 
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6. Relation to conformal field theory 

The symmetric polynomial <E> can be written as a cor- 
relation function of vertex operators V e (z) in a conformal 
field theory (CFT): 35,38 ' 39 

*({*})= lim z^"{V{z O0 )J[V e {z i )) (94) 



V e (which will be called an electron operator) has a form 

V e (z) = i>{z)e^ {z)/ ^ 

where ip is a simple current operator and e 1( t>^ z )/^" is 
the vertex operator on a Gaussian model with a scaling 
dimension h = The scaling dimension h s ^ of ip a (z) = 
[^(z)] a has being calculated from the pattern of zeros 
{S a } in Ref. 11: 

The quasiparticle state $ 7 can also be expressed as a 
correlation function in a CFT: 

$ 1 ({z i }) = lim z 2 S N (Vqiz^V^TJVeizi)) (96) 

Zoo— »00 A A 

( Xol c that $ 7 has a quasiparticle at z = 0.) Here V 1 is 
a quasiparticle operator in CFT which has a form 



Yy(z) =a 7 (z)e i ^ (z)Q - / ^ 7 



(97) 



where a 7 (z) is a "disorder" operator in the CFT gener- 
ated by the simple current operator tp. Different quasi- 
particles labeled by different 7 will correspond to different 
"disorder" operators. 

Let us introduce a quantitative way to characterize the 
quasiparticle operator. We first fuse the quasiparticle 
operator with a electron operators: 

V T , a {z) = V 7 V a = a 7;a (z)e i ^ z )«-/^ 7 

CT 7; a = V-flpa, Q-y:a = Q 1 + a, (98) 

where V a = (V e ) a = </, a e io< W 2 W^. Then, we consider 
the operator product expansion (OPE) of V 7;a with V e 

V e (z)V Tta (w) = (z- w) l °^V Tta+ i(w). (99) 

Let h a , h 7 , and h i;a be the scaling dimensions of V a , Vj, 
and V-y-a respectively. We have 



;Cft _ . _ . _ . 



Since the quasiparticle wave function $ 7 ({zi}) must be a 
single valued function of z^s, this requires that ^ 7; a T must 
be integers. The sequence of integers {^ 7; a T } gives us a 



quantitative way to characterize quasiparticle operators 
V 1 in CFT. 

From the occupation distribution description of the 
quasiparticle 7 introduced in section IIIB, we see that 
a quasiparticle can also be characterized by another 
sequence of integers {^ 7;a }- What is the relation be- 
tween the two sequences of integers, {l^ T } and {l 7;a }, 
that characterize the same set of quasiparticles. From 
eqn. (99), we see that is the order of zeros as we 

move an electron Zi towards a quasiparticle 7 fused with 
a-electrons. Thus is the order of zero -D 7;a ,i in- 

troduced in section III A. From (18) and Si = 0, we 
find that ^ft in the above OPE is given by ^ft = 

S-y,a—S"t,a-i- Thus the two sequences, {^ 7; a T } and {l T , a }, 
are identical {l^ T } — {l T , a }- In the rest of the paper, 
we will drop the superscript CFT in ^ 7; a T - 

Now let us calculate the quasiparticle charge Q 7 (see 
eqn. (97)) from the sequence {l T , a } within the CFT. 
Using l l;a 



■>7, a 



5*7,0-1, we can rewrite (100) as 



*'7;a+l 



S T , a +i — S 7;a + h\. Thus h Ta — 



5 7;a + ahi, where we have used 5 7; o = S* 7 = 0. Using 



^ = hf + ±- = f - ^ , we find 



^7; a ^7 — S^\a ~1~ ^ ^ 

Since ct 7! „ = ay we have 



~2 ~ ~n 



Thus 



(Q 7 + n) 2 -Q n 2 
2v 



a , ( m S n 

o Tn + n\ — 

V 2 n 



m m ^ 



a=l 



which agrees with eqn. (27). 

Let Kf. a be the scaling dimension of a 7 Vv We see that 



n T,a — "7;a 



and 



m 9 (m S-. 



ft 7;a ft 7 -S T , a +^ n ja 



Let 



m 9 /m S-y„ 
we can rewrite the above as 



s -y:a ~ 'J-f.a — + ^ — ja — S 7 ;o+n 



LSC LSC 

"-7;a "-7 



7; a ' 



(101) 



(100) We see that the simple-current part of CFT is determined 



by s 7;a , a = 1, • 



1, only. In particular, if ay and 



07 are related by a simple current operator, ay = a-yV'a, 
then the scaling dimension of ay can be calculated from 
that of 07 : h s 5 — h s ^ + s 7 
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